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PREFACE 


The  research  work  described  in  this  annual  report  was  sponsored 
by  the  Department  of  Defense  under  Project  Themis.  Oklahoma  State 
University  Themis  Weather  Phenomena  Studies  are  administered  under 


contract  No.  DAAB07-68-C-0083. 
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± •  INTRODUCTION 

A.  HISTORY  OF  PROGRAM 

The  College  of  Engineering,  Oklahoma  State  University,  submitted 
Proposal  ER  67-T-21  entitled,  ”A  Center  for  the  Description  of  Environ¬ 
mental  Conditions— Weather  Pnenomena  (Themis  No.  129)''  to  the  Department 
of  Defense  on  May  1,  1967,  in  response  to  the  brochure,  Project  Themis, 
dated  November,  JG56.  Authorization  to  implement  the  program,  t/ith 
limited  funds,  was  given  in  the  form  of  a  letter  contract.  The  effec¬ 
tive  date  of  the  letter  contract  was  October  1,  1967.  The  actual  con¬ 
tract  which  authorized  the  initial  three  years  of  the  program  was 
signed  by  the  contracting  officer  on  January  16,  1968. 

On  January  19,  1968,  a  proposal  for  project  continuation  for  the 
fourth  year  was  submitted  to  the  Department  of  Defense  in  response  to 
the  brochure.  Project  Themis,  dated  November,  1967.  On  June  11,  1968, 
the  contract  between  Oklahoma  State  University  and  the  Department  of 
Defense  was  amended  to  authorize  full  funding  for  the  first  three  years 
of  the  contract.  However,  no  further  notification  of  extension  has 
been  received. 

B.  OBJECTIVES  OF  PROGRAM 

The  proposals  submitted  have  been  based  on  the  following  objectives 

1.  Immediate  objective!  the  determination  of  the  electrot'gnetir 
characteristics  of  severe  weather. 

2.  Intermediate  objective:  the  prediction  of  operationally 
restrictive  and  severe  weather. 

3.  Long-term  objective:  the  modification  of  operationally  re¬ 
strictive  and  severe  weather. 

C.  RA?4 

^nitier  and  definition  of  e;i>ctroemgnetie  signature  assorl- 


atcd  with  specific  severe  weather  occurances  and  the  development  of 
instrumentation  to  detect,  measure,  and  analyze  the  electromagnetic 
energy  present  will  provide  another  technique  for  use  in  short  range 
weather  forecasting.  This  new  method,  along  with  existing  meteoro¬ 
logical  techniques,  will  be  of  benefit  to  the  Department  of  Defense 
in  tactical  operations  of  all  the  services. 

The  modification  of  weather,  if  possible,  will  be  of  benefit  \ 
long-range  military  planning.  The  ability  to  either  predict  or  cause 
a  rainstorm,  for  example,  in  a  particular  area  will  provide  a  tre¬ 
mendous  advantage  in  tactical  operations  planning. 

D.  IMPLEMENTATION  PLANS 

The  program  was  initiated  with  a  three-phased  effort: 

1.  The  development  and  use  of  instrumentation  to  measure  the 
electromagnetic  signature  of  certain  weather  phenomena 
(clear  air  turbulance,  thunderstorms,  tornadoes)  in  airborne 
and  ground-based  laboratories.  (The  correlation  of 
electromagnetic  data  to  existing  meteorological  information 
will  be  made  on  identifiable  weather  phenomena. ) 

2.  The  investigation  of  pattern  recognition  in  the  total  data 
acquisition  program  (sferics,  severe  weather,  existing 
meteorological  information).  (The  mathematical  modeling 
and  information  theory  application  to  both  the  discrete  and 
continuous  data  will  attempt  to  provide  keys  for  the  third 
phase. ) 

3.  The  dete  analysis  and  application  to  identify  indicators 
for  probability  of  specific  phenomena,  determine  possible 
trigger  mechanisms  which  induce  the  phenomena,  end  to  define 
techniques  for  possible  modification-abatement ,  diversion,  or 
intensification  of  the  phenomena . 
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The  development  of  one  or  more  simple  airborne  instrument 
packages  which  can  be  put  in  any  aircraft  for  short  term 
measurement  of  the  severity  of  weairher  immediately  ahead. 
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II.  PROGRAM  MANAGER'S  SUMMARY 

The  program  began  in  October  of  1967  when  funds  became  available, 
equipment  deliveries  began  in  March  of  1968  and  were  completed  in  July 
oi  1968.  Instrumentation  and  shakedown  of  the  research  aircraft  were 
completed  in  August,  1968.  Cooperation  with  the  National  Severe  Storms 
Laboratory  in  the  1966  storm  season  was  not  possible  since  their  data 
gathering  program  was  terminated  on  June  1,  1968,  due  to  the  reduction 
of  their  operational  funds.  Cooperation  with  NSSL  will  be  established 
in  the  spring  of  1969  in  order  to  obtain  close  correlation  between 
meteorological  data  and  electromagnetic  signatures. 

Some  results  with  possible  long-range  meteorological  implications 
have  been  obtained  from  late  summer  storms.  These  results  raise  quest  ions, 
as  to  the  manner  in  which  electrically  charged  clouds  actually  discharge > 
This,  of  course,  is  related  closely  to  the  manner  i.n  which  mass  flow  rf 
air,  both  turbulent  and  non- turbulent ,  occurs  in  a  storm  cloud.  It 
re-emphasises  the  fact  that  future  work  must  now  be  closely  correlated 
with  meteorological  data. 

The  details  of  ths  airborne  instrumentation  package  are  included 
in  a  later  section  of  this  report.  Subject  to  stringent  efforts  to 
tonfirm  some  surprising  results  of  measurements  already  rede,  these 
measurements  indicate  the  following  possibilities: 

1.  At  high  altitudes,  cloud  to  cloud  discharges  often  have 
much  more  horiaontally  polarised  electromagnetic  energy 
than  vertically  polarised  energy.  Cloud  to  ground  discharges, 
without  exception  so  far,  have  massive  energies  in  the 
horisontal  polarisation  a*  well  as  the  vertical  polarixat U  n 
when  measured  at  high  altitudes. 
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2.  The  spectral  energy  distribution  from  flash  to  flash  (of 
the  same  total  discharge)  varies  in  essentially  a  non- 
expected  manner.  This  implies  that  information  on  the 
charge  distribution — and,  therefore,  perhaps  internal  mass 
flow  patterns  and  cloud  structure  information — may  well  be 
obtained  from  horizontally  and  vertically  polarized  spectral 
distribution  information.  The  reason  is  that  polarization 
and  spectral  information  for  a  particular  flash  are  related 
to  the  discharge  path  length  and  path  geometry.  Thus  we 
have  our  first  indication  that  significant  meteorological 
information  (specifically  on  cell  structure)  might  well 
be  obtained  from  the  type  of  measurements  we  are  currently 
making. 

A  detailed  technical  paper,  now  under  preparation,  will  be 
submitted  for  publication  even  though  the  results  are  tentative. 

The  preceding  statements,  though  tentative,  are  rather  far  reaching. 
Therefore  it  is  necessary,  in  the  interests  of  giving  the  reader 
a  somewhat  better  understanding,  to  describe  breifly  the  methods 
of  preliminary  data  reduction  which  have  led  to  these  results.  It 
is  also  necessary  to  describe  three  additional  experiments  which 
will  be  performed  in  the  near  future.  These  experiments  are  designed 
specifically  to  cross  check  these  rather  surprising  results. 

In  preliminary  data  reduction,  an  analog  time  base  expansion 
tachniqua  was  used.  Tigure  1,  shows  the  pkectral  distribution  of  wen 
sive  flashes  of  a  typical  cloud  to  ground  discharge.  Theue  energy 
but  ions  were  obtained  by  taking  horiaontelly  and  vertically  polarized 
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50  Khz  HORIZONTAL 


50  Khz  VERTICAL 

100  Khz  HORIZONTAL 


100  Khz  VERTICAL 


150  Khz  VERTICAL 

200  Khz  HORIZONTAL 

200  Khz  VERTICAL 

250  Khz  HORIZONTAL 

250  Khz  VERTICAL 


FIGURE  1.  CLOUD  TO  GROUND  DISCHARGE 


7 


comb  filter  outputs  which  were  recorded  on  a  multichannel  magnetic  tape 
recorder.  The  tape  was  time  base  expanded  by  a  ratio  of  4:1  and  the 
outputs  displayed  on  a  multichannel  Sanborn  strip  chart.  The  organi¬ 
zation  of  the  signals  on  the  strip  chart  is  such  that  10  KHz  signals  of 
both  horizontal  and  vertical  polarization  are  displayed  adjacent  to  each 
other,  the  same  being  true  for  the  50  KHz,  100  KHz,  150  KHz,  200  KHz, 
and  250  KHz  information.  While  each  individual  flash  is  generally 
identifiable  in  each  record,  the  spectral  distribution  and  polarization 
energies  vary  widely  from  flash  to  flash,  indicating  different  geometries. 
Since  the  visible  part  of  the  discharge  indicated  that  each  flash  has 
the  same  geometry,  one  tentatively  concludes  that  the  intercloud  flashes 
have  radically  different  geometries.  Records  of  this  type  will  be  dis¬ 
cussed  in  considerable  detail  in  forthcoming  papers  to  be  published. 

Preliminary  efforts  are  being  made  to  select  events  of  interest 
from  comb  filter  data  and  then  to  examine  wide  band  records  of  those 
events  in  fine  structure  detail.  Time  base  expansion  in  ratios  of 
1000:1  have  permitted  us  to  examine  both  leader  and  return  of  some 
interesting  events  on  an  analog  basis.  These  results,  though  extremely 
preliminary,  tend  to  confirm  the  previously  indicated  conclusions. 

Within  the  next  few  months,  when  the  digital  reduction  equipment  dis¬ 
cussed  elsewhere  in  this  report  ie  installed,  the  time  base  expended 
wide  band  signals  will  be  fully  digitised.  At  that  print,  complete 
energy  distributions  will  be  determined  through  the  use  of  Fast  Fourier 
transform  techniques  end  our  conclusions  will  be  much  more  firm. 

Severe.!  experiments,  not  previously  anticipated,  will  be  performed 
to  provide  scientific  cross  checks  on  this  matter.  These  ere  necessary 
because  it  ie  not  possible  to  be  absolutely  certain  that  the  signals 
from  each  antenna  have  exactly  the  assumed  polarizations.  At  the  in- 
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ception  of  the  program,  it  was  desired  to  measure  both  polarizations 
to  see  if  there  were  appreciable  differences;  however,  we  did  not 
expect  to  find  the  sharp  differences  encountered,  nor  did  we  expect 
what  differences  we  did  find  to  have  such  significant  meteorological 
implies  Ions.  For  these  reasons,  and  because  the  program  was  running 
late  due  to  delayed  funding,  placement  of  antennas  on  the  aircraft  was 
done  without  static  model  studies.  We  had  considered  reducing  the 
antenna  probes  to  only  one  after  the  first  year.  However,  since  that 
is  now  obviously  out  of  the  question,  a  model  study  experimental  program 
is  being  set  in  motion  to  determine,  as  closely  as  possible,  the  amount 
of  horizontally  polarized  signal  being  received  by  the  vertical  antenna 
and  vice  versa.  Also,  the  aircraft  will  be  flown  to  the  gulf  coast  area 
or  to  the  nearest  available  storm  system  in  the  near  future.  It  will 
then  be  placed  in  a  steep  (60°  or  more)  descending  spiral  from  roughly 
an  altitude  of  20,000  feet  to  10,000  feet  with  all  recorders  running. 

This  violent  maneuver  is  the  only  way  to  alter  the  polarization  char¬ 
acteristics  of  the  vertical  antenna  system  (which  unfortunately  includes 
the  airframe).  It  is  necessary  to  use  a  descending  spiral  rather  than 
a  steep  turn  to  avoid  excessive  G  loads  on  the  aircraft  and  sensitive 
recording  equipment.  It  is  especially  important  to  keep  such  G  loads 
down  since  one  might  encounter  severe  turbulence  at  any  time  near 
r.evoro  storms.  In  the  spring  of  1969,  a  special  high  speed  camera  will 
be  set  up  on  thu  ground  or  in  a  chase  aircraft  to  make  pictures  of  an 
entire  isoleted  storm  structure  at  approximately  250  frames  psr  second. 

At  the  seme  time,  the  research  aircraft  will  be  measuring  tha  storm.  In 
this  way,  we  hope  to  be  able  to  correlate  current  flow  patterns  in  the 
cloud  (by  watching  where  it  glows)  with  individual  cloud  to  ground  flashes, 
while  simulteneoualy  recording  electromagnetic  signals.  Thus  we  hope  to 
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correlate  structure  with  electromagnetic  signatures.  In  these  ways, 
we  believe  we  can  verify  (or  perhaps  deny)  our  findings  to  date. 

In  addition  to  those  experiments  just  discussed,  the  data  gather¬ 
ing  activity  planned  for  1969  is  as  follows.  Beginning  with  the 
1969  storm  season,  data  gathering  flights  will  be  coordinated  with 
the  National  Severe  Storms  Laboratory  program  in  the  Oklahoma  area. 

Thir;  laboratory  makes  minute  by  minute  photographic  records  of  contour 
radar  pictures  covering  about  a  200  mile  radius  of  Stillwater.  This 
information  gives  good  data  on  moisture  content  and  severe  storm  renters. 
NSSL  also  gathers  considerable  amounts  of  other  meteorological  data 
in  their  flight  and  ground  record  programs.  There  is  no  better  way 
to  tie  our  electromagnetic  measurements  to  pertinent  and  accurate 
meteorological  data;  thus,  from  about  February  1  to  June  1,  our  entire 
activity  will  be  in  the  Oklahoma  area. 

Tor  perhaps  a  total  of  ten  days  in  July  or  August  of  1969,  we 
plan  to  make  measurements  at  the  Organ  Peak  Laboratory  at  White  Sands. 

The  opportunities  for  getting  data  of  repeatable  significance  is 
perhaps  greater  here  than  anywhere  else.  Our  aircraft  will  be  available 
on  Condron  field  for  immediate  flight  when  a  storm  cloud  develops 
over  Organ  Peak.  Storms  occur  on  roughly  thirty  percent  of  the  days 
during  July  and  August.  Thus,  we  should  be  able  to  get  meaningful 
correlation  between  electromagnetic  signatures  and  the  size  and  state 
of  maturity  of  an  isolated  storm  for  many  essentially  similar  storms. 

We  will  simultaneously  have  observers  on  Organ  Peak  in  communicat ion 
with  the  aircraft  so  that  full  meteorological  information  is  made  part 
of  tho  minute  by  minute  petwanent  record. 

One  piece  of  equipment,  a  field  mill,  will  be  added  to  the  aircra** 
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:;o  that  electrostat r  field  becomes  a  part  of  the  permanent  record. 

This  addition  is  also  believed  to  be  necessary  because  the  aircraft 
has  encountered  at  least  o.  e  condition  of  extreme  electrostatic  fields 
that:  disabled  much  of  the  other  equipment.  This  situation  is  discussed 
elsewhere  in  this  report. 

Also,  in  another  section  of  this  report,  a  discussion  is  given 
of  the  digitizing  equipment  that  will  be  used  primarily  on  the  wide 
band  analog  data.  With  this  capability  in  full  operation,  it  will  b? 
possible  to  examine  in  much  greater  detail  tne  fine  structure  of  the 
electromagnetic  signature  using  digital  Fast  Fourier  transform  tech¬ 
niques.  The  importance  of  this  type  of  analysis  is  becoming  more 
obvious  with  the  results  already  obtained. 

The  pattern  recognition  group  has  been  greatly  hindered  by  the 
l act  that  data  gathering  has  been  so  far  behind  schedule.  Their  re¬ 
port  is  detailed  in  another  section  of  this  document.  They  have  been 
making  headway  with  the  data  that  is  now  available,  and  have  already 
demonstrated  the  digital  techniques  they  intend  to  employ  by  using 
limited  digitizing  equipment  currently  available. 

The  cloud  modeling  group  report  is  also  included.  This  report 
essentially  terminates  that  phase  of  the  program.  Our  major  effo:  s 
I rom  here  on  are  to  be  concentrated  in  the  data  gathering,  reduction, 
and  pattern  recognition  areas.  This  decision  to  shift  emphasis  is 
dictated  by  r.erious  need  to  very  closely  correlate  electromagnetic 
signatures  with  minute  by  minute  meteorological  data. 

In  summary,  ths  achievements  to  date  are  listed  ss  follow;* 

1.  Design,  fsbrication,  installation,  shakedown,  and  operation 
of  a  sophisticated  airborne  instrumentation  package. 

2.  Development  and  checkout  of  analog  data  reduction  methods 
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whose  application  to  somewhat  meager  data  has  led  to  interesting 
tentative  conclusions  and  has  called  for  several  new  experiments. 

3.  Development  and  checkout  of  digital  reduction  methods  that, 
when  combined  with  the  use  of  equipment  soon  to  be  installed, 
will  allow  analysis  of  already  discussed  new  phenomena  in 
much  more  highly  resolved  detail. 

4.  Development  and  initial  application  of  pattern  recognition 
methods  directed  towards  correlation  of  meteorological  data 
with  electromagnetic  signatures,  the  ultimate  objective  of 
the  program. 

In  the  course  of  carrying  out  this  work  for  the  past  year,  many 
conferences  have  been  held  with  other  researchers  in  the  field.  Some 
of  these  are  listed  below. 


December  18-19,  1967  Location:  White  Sands  Proving  Grounds 

New  Mexico 

Subject:  Instrumentation  of  Research 

Aircraft 

Participants:  Marvin  Diamond  -  White  Sands 
E.  A.  Blomerth  -  White  Sands 
William  L.  Hughe*  -  OSU 
Arthur  H.  Breipohl  -  OSU 

March  8,  1968  Location:  National  Center  for  Atmospheric 

Research,  Boulder,  Colorado 

Subject:  Electromagnetic  Radiation 

Participants:  Dwayne  Sarrtor  -  NCAR 

William  L.  Hughes  -  OSU 
Arthur  M.  Breipohl  -  OSU 
Victor  W.  Bolie  -  OSU 

dune  7,  1968  Location:  National  Severe  Storm 

Laboratory,  Norman,  Oklahoma 

Subjrct:  Coordination  of  NS  SI.  and  OSU 

Weather  Research  Project* 

Participants:  Edwin  Kessler  -  NSSL 
Gilbert  Kinter  -  NSSL 


1. 

{ 

i 

.  - 

-i 

12 

Gene  Lee  -  NSSL 

William  L.  Hughes  -  OSU 

Arthur  M.  Breipohl  -  OSU 

•July  1,  1968 

Location: 

Oklahoma  State  University 

Stillwater,  Oklahoma 

Subject: 

Review  Themis  Weather 

Phenomena  Program 

Participants: 

G/lbert  Kinzer  -  NSSL 

William  L.  Hughes  -  OSU 

Arthur  M.  Breipohl  -  OSU 

Victor  W.  Bolie  -  OSU 

July  2,  1968 

Location: 

Oklahoma  State  University 

Stillwater,  Oklahoma 

Subject: 

Review  Themis  Weather 

Phenomena  Program 

Participants: 

W.  L.  Taylor  -  ESSA 

William  L.  Hughes  -  OSU 

Arthur  M.  Breipohl  -  OSU 

Victor  W.  Bolie  -  OSU 

July  10,  1968 

Location: 

Oklahoma  State  University 

Stillwater,  Oklahoma 

Subject: 

Aircraft  Instrumentation  and 

Radar  Equipment 

Participants: 

Col.  Robert  W.  Vincent, 

Air  Weather  Service 

Robert  L.  Overton  -  OSU 

July  14-15,  1968 

Location: 

Oklahoma  State  University 

Subject: 

Site  Visit  -  Department  of 

Defense  Monitors 

Participants: 

Themis  Weather  Phenomena 

Project  Staff 

Trances  L.  WheJon  -  P.O.D. 

Wilbur  Paulson  -  A.T.  Cambridge 

E.  A.  Blomerth  -  White  Sands 

R.  Loveland  -  White  Sands 

August  5-9,  1968 

location: 

White  Sands  Missile  Range, 

New  Mexico 

1 

Cub jecv : 

Electromagnetic  Radiation  Data 

Gathering  flights  in  Single  Cell 

Stormn  and  Meetings  with  Atmo¬ 
spheric  Science  Laboratory  Personnel. 

1 

1 _ 

1? 


September  17,  19o8 


November  13,  1968 


Participants;  William  L.  Hughes  -  OSU 
Gerald  Stotts  -  OSU 
Keith  L.  Mackey  -  OSU 
Russell  L.  Babb  -  OSU 
Marvin  Diamond  -  White  Sands 
Alex  Blomerth  -  White  Sands 


Location:  Institute*  of  Atmospheric 

Physics 

Subject:  Discuss  Themis  Project 

Obtain  Consultant  for  Project 


Participants:  William  L.  Hughes  -  OSU 

L.  J.  Battan  -  Institute  of 
Atmospheric  Physics 

Location:  Oklahoma  State  University 

Stillwater,  Oklahoma 


Subject:  Instrumentation  of  Research 

Aircraft 


Participants:  William  L.  Hughes  -  OSU 

George  Kirkpatrick  -  Uni¬ 
versity  of  Syracuse 


A  serious  deficiency  in  the  program  has  been  the  lack  of  on-the- 
spot  meteorological  expertise.  We  have  hired  Dr.  Louis  Battan  as  a 
consultant  to  the  project  and  are  making  major  efforts  towards  hiring 
a  full  time  meteorologist.  One  offer  was  made  to  a  highly  qualified 
Ph.D.,  but  after  several  weeks  of  consideration,  he  has  declined.  We 
are  vigorously  pursuing  several  other  leads  at  the  moment.  In  the 
meantime.  Dr.  Battan  has  exhibited  considerable  interest  in  the  project 
and  has  made  many  helpful  suggestions  about  the  program. 

The  current  work,  during  the  lull  in  available  storms,  has  been 
directed  toward  documenting  what  has  already  been  done. 

One  Ph.D.  thesis  has  resulted  from  the  program: 

1.  An  Analysis  of  Charge  Transport  and  Equilibrium  in  the 
Global  Atmosphere. 


by:  Edward  L.  Shreve,  School  of  Electrical  Engineering. 
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Three  technical  reports  have  been  written: 

2.  a.  Preliminary  Survey  of  the  Present  State  of 

Knowledge  of  Atmospheric  Electricity. 

by:  Edward  L.  Shreve,  School  of  Electrical  Engineering. 

b.  Preliminary  Survey  of  the  Present  State  of  Knowledge 

on  Lightning  ” 

by:  Eddy  J.  Milanes,  School  of  Electrical  Engineering. 

c.  Preliminary  Survey  on  the  Present  State  o 2  Knowledge  of 
the  Spectra  Prom  Lightning-Produced  Plasma 

by:  William  G.  Robinson,  Departir vnt  of  Physics 
One  technical  paper  has  been  submitted  for  publication: 

3.  Theoretical  Derivation  of  Atmospheric  Electrical  Properties 

by:  Edward  L.  Shreve,  School  of  Electrical  Engineering 
submitted  to:  Journal  of  Geophysical  Research 
Four  technical  papers,  are  being  written  for  submission  for 
put)  li  cat  ion.  Tentative  titles  are: 

4.  a.  "Electric  Field  and  Potential  Around  Some  Thunderstorm 

Models" 

b.  "Electromagnetic  Radiation  Related  to  Current  in  Lightning 
Discharges" 

c.  "Differences  Between  Horizontally  and  Vertically  Polarir.ed 
Sferic  Signals  at  High  Altitudes" 

d.  "Analog  Method  for  Preliminary  Reduction  of  Sferic  Signals" 
Copies  of  the  thesis,  three  technical  reports,  and  the  technical  pat .  r 


sulxnitied  for  publication  are  included  in  Appendix  B  (a  separately  K-und 
addendum  to  thin  report). 
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III.  DAI A  GATHERING 

This  section  traces  the  development  of  the  airborne  data  gathering 
assembly  throughout  the  reporting  period.  Included  are  the  results 
of  system  performance  tests  and  the  self-checks  using  recorded  weather 
data.  In  addition,  the  airborne  data  taken  during  the  past  year  is 
discussed  in  the  context  of  system  evaluation. 

The  test  bed  for  the  measurement  apparatus  is  a  D-18S  Reechcraft 
airplane.  Because  the  singular  mission  of  the  aircraft  is  data  gather¬ 
ing  of  weather  informatior ,  approval  was  obtained  from  the  Federal 
Aviation  Administration  tc  operate  in  tne  "Restricted"  category.  Since 
the  aircraft  was  procured  in  a  used  condition,  a  renovation  was  initiated. 
Concurrent  with  the  renovation  of  the  aircraft,  provisions  were  made 
for  t.he  mounting  of  a  vertical  whip  antenna  atop  the  fuselage  and  a 
tow  type  antenna  extending  horizontally  from  the  tail  of  the  aircraft. 

A  take-up  or  release  reel  has  been  incorporated  in  the  horizontal  an¬ 
tenna  to  facilitate  the  take  off  and  landing  of  the  aircraft.  The 
aircraft  appears  as  shown  in  Figure  2.  The  measurement  assembly  is 
mounted  in  the  airplane  as  shown  in  Figure  3.  The  rack  consists  ot  the 
following  commercial  equipment: 

1.  2  Sangemo  multichannel  recorders,  1  fourteen-track  FH  and  1 
five-track  wide  band. 

2.  Hewlett-Packard  oscilloscope. 

3.  Special  Products  HWV  receiver. 

«*.  Topes  power  inverter. 

5.  American  Radio  Corporation  transceiver. 

6.  Military  UHf  transceiver. 

Additional  equipment  developed  expressly  for  this  project  ar*: 

1.  2  comb  filter  penels. 


2. 


18 


Remote  program  panel  for  Sangamo  14  channel  recorder. 

3.  DC  power  supply  filter  panel. 

4.  Monitor  panel. 

The  power  distribution  schematic  for  the  assembly  is  shown  in 
Figure  4.  Ail  DC  power  is  derived  directly  from  the  aircraft  generation 
system.  AC  power  requirements  of  the  oscilloscope  and  the  remote 
program  panel  are  met  by  the  use  of  the  Topaz  power  inverter.  In  the 
early  stages  of  fabrication  of  the  data  gathering  assembly,  tests 
revealed  that  interference  from  the  inverter  was  being  coupled  into 
the  aircraft  DC  generation  system.  This  effect  necessitated  the 
incorporation  of  specially  designed  filters.  Through  the  use  of  these 
filters,  the  interference  has  been  sufficiently  reduced. 

The  measurement  channels  are  fed  by  the  two  antennas  referred  to 
earlier.  This  allows  the  simultaneous  recording  of  vertically  and 
horizontally  polarized  signals.  The  vertical  signal  flow  path  is  shown 
in  Figure  b.  The  signal  gathered  at  the  antenna  is  fed  into  a  low 
par'  filter  and  then  through  a  pro amp/cathode  follower  circuit  (see 
Figures  b  and  /).  Nuvistors  ire  being  used  in  this  stage  duo  to  the 
occurrence  of  largo  transients  and  static  charges  which  preclude  the 
use  of  solid  state  devices.  The  low  pass  filter  became  necessary  when 
it  w.v-.  discovered  that  grid  detection  of  the  ship's  air  to  ground 
•  ••  •mmun  ieat  ion  signal  war.  occurring,  at  the  fii  t  cathode  follower  ■  t  , 

I  be  detected  audio  signal  blanked  tin;  ins  t  rumen  tat  ion  during  radio 

I I  itiam i ss ion .  because  the  carrier  of  the  air-ground  common i cut  ion 
th.intiel  war;  well  beyond  the  cut-off  of  the  measure-nont  system,  it  w.i. 
go-  tide  to  insert  the  low  par;  filter  without  appreciably  affect  iug. 

tin-  amplitude  «.r  flatnc.o  of  'he  do  iron  sp«v tr>jn.  Fu;-  \...  .i  -  ignis'  *  -mf 
bit'  jn  -  jgfial  strength  id.  ;i  vert  i » •  1 1  u\  l  :ie»  .  ■  »i  •  h  »r»>--  • , 
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r ! G U R E  6  VERTICAL  PREAMP/ CATHODE  FOLLOWER 


it  also  became  necessary  to  include  a  preamp  stage  in  the  vertical  chan 
nei.  This  stage  has  a  pain  of  about  10  lbs  and  raises  the  vertically 
polarized  signal  to  the  comparable  range  of  the  horizontal  signal. 
Included  in  the  preamp/cat!  ode  follower  circuit  is  a  two-stage  cathode 
follower  network  connect •'  in  cascade.  The  first  stage  acts  as  a 
buffer  to  absorb  high  amplitude  transients  and  static  changes  appearing 
on  the  antenna.  Owing  to  the  quiescent  bias  on  the  second  stage,  sat¬ 
uration  of  the  first  stage  is  not  fed  through  the  system.  In  addition, 
the  recovery  time  of  the  system  is  considerably  reduced. 

From  the  preamp/cathode  follower  stage  the  signals  are  fed  into 
the  comb  filter  panel.  Each  comb  filter  panel  consists  of  six  bandpass 
detector  channels  and  one  broad  band  channel.  The  circuit  diagrams 
for  typical  cards  of  each  type  are  shown  in  rigurer,  8  and  0.  Center- 
frequencies  of  the  bandpass  cards  are  10  KHz,  50  KHz,  100  KHz,  150  KHz, 
200  KHz,  and  250  KHz.  The  pass  band  is  500  Hz  each  side  of  center- 
frequency.  The  filtering  is  accomplished  through  the  use  of  opera! Lon a 
amplifiers  which  results  in  a  gain  cf  100  over  the  pass  band  The 
dynamic  amplitude  range  of  the  comb  filters  is  0-1  millivolt  peak 
input  and  0-100  peak  millivolt  output.  Each  filter  module  performs  an 
additional  function,  namely,  envelope  detection  with  a  detection  time 
constant  of  0.001  seconds.  Thus,  the  signal  appearing  at  the  output 
of  each  channel  in  the  envelope  of  the  excersions  about  the  center 
frequency.  The  broad  band  channel  does  not  have  a  detector  system. 
Father,  the  entire  spectrum  from  roughly  3  KHz  to  300  KHz  is  ampliti.  , 
by  a  factor  of  30  and  then  recorded  with  no  further  processing. 

The  output  of  the  comb  filter  panels  is  connected  through  the 
monitor  panel  to  the  recorders.  The  monitor  panel  is  essentially  a 
feedthrough  device  where  the  nignal  paths  are  selectively  tapped  by  t 
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FIGURE  8  NARROW  BAND  AMPLIFIER 
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FIGURE  9  TYPICAL  BROADBAND  AMPLIFIER  -  VERTICAL  AND  HORIZONTAL 


23  position  rotary  switch. 

The  switch  wiper  is  connected  to  the  oscilloscope.  In  addition 
to  the  comb  filter  outputs,  a  tap  has  been  provided  for  the  monitoring 
ol  each  jf  the  comb  filter  input  signals.  Thus,  a  real  time  display 
of  the  d  ita  before  and  after  processing  is  continually  available. 

This  allows  the  prediction  of  optimal  recording  levels. 

The  detected  bandpass  data  is  recorded  on  12  channels  of  one  of 
the  14  channel  Cangamo  recorders.  This  recorder  is  equipped  with  FM 
recording  modules  to  accommodate  the  low  frequency  characteristics  of 
the  detected  uata.  During  the  data  gathering  mission,  this  recorder 
operates  continually  at  a  tape  speed  of  7  1/2  inches  per  second.  The 
broad  band  data  output  of  the  comb  filter  panel  is  fed  into  4  channels 
of  a  second  Sangatno  recorder.  Each  of  the  broad  band  outputs  is  re¬ 
corded  on  an  FM  module  and  a  direct  module  simultaneously.  This  re¬ 
corder  is  operated  for  1  minute  every  IS  minutes  at  a  tape  speed  of  60 
inches  per  second.  The  necessity  for  the  parallel  recording  of  each 
channel  on  two  tracks  results  from  the  fact  that  the  data  will  be 
recovered  at  1  7/8  inches  per  second  to  gain  a  time  base  expansion. 

In  this  process,  the  lower  frequency  spectrum  of  the  direct  record 
module  will  be  lost  due  to  the  frequency  response  of  the  reproducing 
head.  However,  this  data  will  be  recovered  from  the  FM  replay  module. 

In  addition  to  the  weather  data,  both  recorders  are  recording  a  time 
reference  trom  the  WWV  receiver.  This  time  reference  is  applied  to 
voice  edgetracks  in  the  recorders.  The  recorder  storing  the  detected 
outputs  has  a  mc-dification  introduced  in  the  edgatrack  which  atso  allows 
voice  cueing.  In  addition,  provisions  have  been  mad#  to  record  data 
f»o*  a  ground  station  onto  a  second  voice  track  contained  in  thin  re 
confer.  The  output  of  a  VHP  transceiver  will  be  fed  to  this  edgrtra  !>  , 
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The  intent  is  to  correlate  the  storm's  electrical  data  with  meteorological 
data  observed  from  the  ground  or  from  another  chase  aircraft. 

Because  of  the  inconvenience  experienced  by  the  airborne  tech¬ 
nician:;  in  the  task  of  intermittently  operating  the  wideband  recorder, 
a  iemote  Program  Panel  was  developed.  This  panel  consists  of  an 
industrial  timer  with  a  cycle  of  15  minutes  and  a  variable  period  of 
actuation.  Associated  digital  logic  translates  the  activate  period 
into  "record"  and  "stop"  commands  to  the  recorder. 

An  additional  transceiver  installed  in  the  aircraft  is  a  UHf 
communicator.  This  was  required  by  range  safetv  personnel  while  the 
airplane  was  being  operated  in  the  vicinity  of  White  Sands. 

Following  the  fabrication  of  the  measurement  system  and  static 
bench  tests,  the  assembly  was  sub-jected  to  airborne  noise  tests.  The 
first  three  flights  were  primarily  for  the  purpose  of  shaking  down 
the  system.  The  results  of  these  flights  resulted  in  several  minor 
modifications.  For  example,  the  horizontal  antenna,  which  is  a 
trailing  wire  antenna,  was  found  to  be  flapping  excessiv  ly  during 
flight.  This  resulted  in  considerable  perturbations  in  the  horizontal 
channel.  This  problem  was  corrected  by  increasing  the  weight  of  the 
nock  at  the  end  of  the  antenna.  The  vertical  antenna  was  found  to  hr 
picking  up  signals  from  the  a:icraft's  air-to-ground  communications 
equipment.  The  antenna  was  moved  forward  several  feet  to  increase  it:, 
distance  from  the  communications  antennas.  A  filter  w us  added  to  th.> 
system  to  filter  out  the  redi<>  frequency  signals.  Uter  it  was  found 
the  preamplifier  in  the  horizontal  channel  was  not  required.  A  pro* -dm.- 
wa?  determined  for  the  equalization  of  the  gains  in  the  two  channels. 

This  is  done  by  equalizing  the  white  noise  levels  in  the  wideband 
channels  by  adjusting  the  length  of  the  horizontal  antenna  while  in  «iir*. 
(The  horizontal  antenna  may  be  reeled  in  or  out  bv  the  technician  ihtr! 


the  aircraft  merely  by  operating  a  switch  on  the  data  gathering  panel.) 

nights  4  and  l>  were  undertaken  after  the  system  was  completely 
operational,  blight  4  was  very  significant  because  it  was  the  first 
(light  iti  which  useful  data  was  obtained.  The  storm  observed  was  a 
reasonably  isolated  cell  of  moderate  intensity  and  of  moderate  physical 
si/.e.  The  aircraft  remained  in  clear  air  at  all  times  during  the  flight, 
flight  5  differed  from  flight  4  in  several  aspects.  The  cells  in  this 
instance  were  not  as  well  isolated  as  the  previous  flight.  The  cells 
were  quite  large,  extending  perhaps  to  fifty  or  sixty  thousand  feet. 

A  large  shelf  extended  many  miles  out  from  the  cells.  Much  of  the 
flying  was  done  under  this  shelf.  One  interesting  phenomena  was  ob¬ 
served  while  flying  under  this  shelf.  Large  static  charges  accumulated 
on  both  the  horizontal  and  vertical  antennas.  These  charge  accumu- 
1 <t ions  and  the  subsequent  discharges  caused  the  comb  filter  inputs 
to  saturate  several  times  during  the  flight. 

Because  the  storm  season  was  virtually  expended  in  the  Oklahoma 
area,  nights  7,  8,  and  9  were  undertaken  in  the  White  Sands  area.  The 
data  gathered  during  these  flights  are  particularly  valuable  since  the 
storm  cells  were  well  isolated  and  of  moderate  to  bigh  electrical  in¬ 
tensity.  In  flights  8  and  9,  the  visual  characteristics  of  the  elec¬ 
trical  activity  were  entered  on  the  edgetrack.  The  reduction  of  this 
data  Is  currently  in  progress.  One  particularly  ancouraging  result 
of  this  reduction  was  the  outcome  of  a  uelt-check.  The  broad  band 
data  was  replayed  through  a  comb  filter  on  test  bench  power.  The 
output  of  the  filter  was  an  exact  match  with  the  narrow  bant  data  ;•»>. 
forded  earlier  in  the  data  gathering.  Hence,  the  repea? ibi Uty  of  the 
system  holds. 

In  conclusion,  it  i*.  selt  that  the  data  gathering  a-r-eoblv  i-;  .>m- 

ple.ely  operational  and  further  *m  jor  modif  i<‘*l  ions  t'>  the  existing 


/  ’  1 

functional  arraignment  arc  rot  contemplated.  The  only  change  anticipat'd 
i  to  install  a  field  mill  so  that  static  fields  can  If  recorded. 
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IV.  PATTERN  RECOGNITION 

Th?  purpose  of  the  pattern  recognition  research  is  to  detennine 
il  the  oi.ectronwip.net ic  signature  of  a  storm  contains  information  which 
can  be  used  to  identify  certain  meteorological  phenomenon.  Pattern 
recognition  may  be  accomplished  by  simply  observing  the  signals  re¬ 
corded  by  the  measuring  instruments;  however,  this  phase  of  the  re¬ 
search  is  based  on  the  assumption  that  there  may  be  patterns  in  the 
received  signal  which  can  be  detected  only  by  automatic  (computer) 
pattern  recognition.  With  this  premise,  all  pattern  recognition 
efforts  are  based  upon  auU matically  reducing  the  recorded  signal  to 
a  usable  form  for  the  compute. 

Research  to  the  present  time  has  been  divided  into  two  cate¬ 
gories,  parametric  and  nonparametric  pattc  -n  recognition.  In  this 
pattern  recognition  problem,  the  major  problem  is  to  reduce  the  number 
of  measurements  to  a  size  that  can  be  used  bv  the  computer  and  yet 
retain  the  significant  information.  In  the  parametric  pattern  rec¬ 
ognition,  we  hope  to  select  measurements  that  seem  reasonable  and 
let  the  computer  determine  if  the  selected  measurements  are  sufficient 
to  predict  meteorological  phenomenon. 

The  following  section'-,  describe  these  two  phases  of  this  research. 
A,  PARAMETRIC  PATTERN  RECOGNITION 

In  thin  phase  of  the  project,  a  model,  based  on  the  sferic 
measurements,  is  being  derived  to  Hopefully  aid  in  the  classification 
of  a  given  storm.  There  is  no  intention  to  model  exactly  the  electro¬ 
magnetic  phenomena  occurring  within  the  storm.  What  is  desired  is  a 
model  that  will  represent  th«  average  effect  of  these  phenomena  in 
'•u<  h  a  manner  that  it  may  te  used  successfully  in  a  pattern  classi¬ 
fication  process.  The  procedure  being  followed  t  •  derive  this  modei 
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may  be  summarized  us  follows: 

1.  Obtain  a  baric  model  for  a  general  sferic  generating 
discharge. 

2.  Relate  this  model  to  observation:;  at  the  antenna. 

3.  Use  the  sferic  measurements  at  the  antenna  to  improve  the 
model. 

4.  Ultimately  use  the  refined  model  to  classify  storms. 

At  this  point,  a  basic  model  has  beer  assumed  that  consists  of 
an  electric  and  a  magnetic  dipole.  These  dipoles  have  been  related 
to  the  sferics  measured  at  the  antenna  through  conventional  electro¬ 
magnetic  field  calculations.  Sferic  recordings  obtained  on  past  storms 
are  now  being  used  to  determine  the  form  of  the  improved  model.  The 
results  obtained  under  each  of  the  outlined  steps  in  the  procedure 
will  be  discussed  in  the  following  four  sections,  respectively. 

Basic  Model 

An  exact  model  of  the  individual  electrical  activities  occurring 
within  the  storm  will  not  be  attempted.  Tor  example,  lightning  strokes 
obviously  contribute  to  the  measured  sferics;  however,  lightning 
strokes  will  not  be  modeled  as  such.  One  reason  for  this  approach 
is  the  difficulty  in  modeling  due  to  the  randomness  associated  with 
the  paths,  times  of  occurrence,  waveforms,  and  branching  conditions 
of  *he  strokes.  However,  most  important  is  the  possibility  that  the 
formulation  of  a  modal  in  terms  of  lightning  strokes  and  other  familiar 
discharges  may  neglect  3ome  type  of  discharge  that  is  not  familiar  bur 
contributes  to  the  sferic  measurements  and  contain?  information  in¬ 
valuable  for  storm  classification.  Thus,  rather  than  requiring  the 
basic  model  to  represent  several  particular  electrical  phenomena  within 
the  storm,  it  is  required  that  the  basic  model  represent  a  general 
sferic  generating  element.  An  electric  dipole  and  a  magnetic  dipole 
have  been  chosen  to  represent  this  element.  Later  in  the  procedure. 


several  of  these  dipoles  will  be  combined  into  an  improved  model  such 
that  the  sferics  generated  by  this  improved  model  will  agree  with  the 
ul erics  measured. 

The  assumption  of  dipoles  for  the  basic  model  may  seem  somewhat 
over  simplified.  However,  as  stated  earlier,  it  is  not  desired  to 
obtain  an  exact  model  of  the  actual  discharges  which  are  occurring. 

Tt  is  intended  to  formulate  a  model  which  represents,  on  the  average, 
the  sferics  as  measured  by  the  antenna  so  that  a  pattern  classification 
process  may  be  applied.  Since  the  mechanisms  which  generate  a  par¬ 
ticular  sferic  are  not  unique,  the  model  may  be  as  flexible  as  desired 
us  long  as  the  sferics  it  generates  correspond  to  the  measurements 
at  the  antenna. 

The  remainder  of  this  section  will  be  devoted  to  describing  the 
form  in  which  the  electric  dipole  is  to  be  used  in  the  basic  model. 

The  notation  for  the  magnetic  dipole  will  not  be  discussed  other  than 
noting  that  it  has  a  form  identical  to  the  electric  dipole  with  the 
subscripts  changed  from  e  to  m  on  all  variables  and  parameters. 

In  reference  to  figure  10,  the  dipole  is  assumed  to  be  located 
at:  the  random  point  iX  ,  Y  Z  )  with  orientation  (1  ,  m  ,  n  ).  The 
<>l> MTVor  is  loi.ited  at  the  point  (x,  y,  r.).  Since  the  rferies  are  not 
expected  to  occur  at  fixed  intervals  of  time,  neither  should  radiation 
fr<»m  t tie  dipole  occur  in  this  manner.  Thus,  it  is  assumed  that  the 
dipole  will  begin  radiation  at  the  random  times  with  current 
waveforms  The  current  expression  for  the  dipole  becomes 
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This  process 


is  commonly  known  as  "shot  noise". 
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Figure  10.  COORDINATE  SYSTEM 

The  basic  model  explicitly  involves  the  orientation  parameters 

1  ,  m  ,  n  which  must  be  learned.  In  addition  to  these,  the  form  of 
e  e  e 

the  distribution  functions  must  be  derived  for  the  assumed  random 
point  of  occurrence  in  space  and  time.  These  distribution  functions 
will  undoubtedly  contain  additional  parameters  which  must  be  learned. 
Also,  it  is  assumed  that  the  waveform  ;»j.(t)  is  time  harmonic  and  mav 
he  expressed  in  a  Tourier  series  which  introduces  more  parameters  in 
terms  of  the  Fourier  coefficients. 

The  learning  of  these  distributions  along  with  all  parameters 
introduced  will  he  accomplished  through  the  sferic  measurements  it  ni 
antenna.  This  :s  the  topic  for  discussion  in  section  t-C. 

t!efor«  the  antenna  data  is  used  to  learn  the  distributions  and 
p»r.»««iteru  of  the  basic  model  and  to  expand  this  model,  the  radiation 
from  the  basic  model  must  be  < alculated  so  that  it  is  related  to  the 
antenna.  These  calculations  will  he  discussed  in  the  following  so. "  i 
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Field  Calculation?, 

The  electric  dipole  may  bo  represented  by  a  positive  time  varying 

charge,  Q(t),  separated  from  an  equal  but  opposite  polarity  charge 

-Q(t),  by  -  distance  AL  along  the  unit  vector  (1  ,  m  ,  n  ).  (See 

e  e  e 

Figure  11.)  Mathematically,  the  dipole  may  be  expressed  by  the  po¬ 
larization  vector 

~P(t)  =  Q(t)  ALT  , 

where  a  is  the  unit  vector  (1  ,  m  ,  n  ). 
e  e’  e’  e 


i-  Q(t) 

AL  / 

/  / 

°  -0(t) 


Figure  11.  DIPOLE  REPRESENTATION 


The  field  distribution  will  be  calculated  by  determining  the 
He’' :  z  vector  at  the  observation  point  (x,  y,  z).  The  Hertz  vector 
resulting  from  the  above  polarization  vector  located  at  the  point 
(x' ,  y',  z*  )  is  given  by 


Il(x,  y,  z,  t)  = 


-  r/v)  —  dv' 
r 


2  2  21/2 

where  r  is  the  distance  ({x  -  x')  ♦  (y  -  y * )  ♦  (z  -  z*)  j  and 
v  is  the  propagation  velocity  of  the  wave.  In  these  calculations,  AL 
will  be  chosen  such  that  the  Hertz  vector  may  be  approximated  as 

!l<x,  y,  z,  t)  i  t  -  r/v)  , 

O 
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or 

H(x,  y,  2,  t)  =  r  Q(t  -  r/v)  AL“ae  . 

o 

The  field  conditions  are  given  .in  terms  of  the  Hertz  vector  by 
the  equations 

“?=  V(V  n)  -  ue  , 

at 


<iml 


B 


ye 


Assuming  the  dipole  is  time  harmonic,  the  above  two  equations  beet  me 


E  *  V(V  it)  +  uew  n 


and 

B  =  juiye  7XH  . 

Before  continuing  with  the  procedure,  two  reasons  for  using  the 
Hertz  vector  approach  to  the  field  calculations  will  be  noted,  first, 
the  dipole  current  distribution  is  given  by  J  =  3P/3t.  Thus,  rather 
than  requiring  the  current  distribution  to  be  known,  the  chosen  formu¬ 
lation  needs  only  its  time  average.  Second,  the  field  conditions  fot 
the  magnetic  dipole  may  be  obtained  with  ease  from  the  field  conditions 
of  the  electric  dipolo.  This  will  be  discussed  at  the  end  of  the  section. 

In  determining  the  field  conditions  from  the  equations  for  f  tnd  »» 
the  calculations  will  not  he  limited  to  the  far  field.  Kraus  1111  m  it. 

^References  cited  are  shown  in  Appendix  A. 


that  the  boundary  between  the  near  and  fe.r  fields  for  a  dipole  antenna 

2 

may  bo  taken  to  be  at  a  radius  R  =  2b  /X,  where  L  is  the  length  of  the 
antenna  ana  X  ir.  tne  wavelength  of  the  radiation  frequency.  If  the 
model  were  required  to  represent  a  cloud  to  ground  stroke  one  mile  in 
h<  irlit,  a  sample  R  may  he  calculated  using  the  frequencies  of  10  KH. 
and  2 SO  KH. ,  which  are  the  extremes  of  the  narrow  band  recordings. 
Tiio.oe  frequencies  yield  R  =  0.41  miles  for  10  KH.  and  R  =  10.7  miles 
tor  250  KH.  Thus,  at  the  higher  frequencies  we  are  not  guaranteed 
t  ir  field  conditions,  and  the  calculations  must  include  all  components 
of  the  field. 

To  account  for  the  boundary  conditions  at  the  earth's  surface, 
imago  theory  is  assumed  to  hold.  Under  this  assumption,  the  effect 
of  the  earth's  surface  may  be  represented  by  the  original  dipole's 
image.  This  consists  of  a  dipole  identical  to  the  original  dipole  and 
located  a  distance  of  2h  directly  below  the  original  dipole  with 
orientation  (-1(1,  -me,  nc>.  Here  it  is  assumed  that  the  original 
dipole  was  located  at  a  height  h  above  the  earth's  surface.  The  wave 
guide  effect  produced  by  the  earth's  surface  and  the  ionosphere,  and 
the  ground  wave  effect  have  briefly  been  investigated  and  believed  to 
be  negligible.  Thus,  the  problem  of  calculating  the  radiation  from 
■in  electric  dipole  located  a  height  h  above  the  earth's  surface  will 
be  approximated  by  calculating  the  radiation  from  the  original  dipole 
and  its  image  in  free  space. 

for  the  pjesent  calculation,  the  dipole  will  be  located  at  the 
point  (o,  o,  h)  and  its  image  at  (o,  o,  -h);  the  observer  is  at  (x,  y. 
The  Hortz  vector  a?  the  ob'-'^vation  point  is  given  by 


tn«  ^  ^  0.  Q,  Q. 

(-7T2)  n  «  a  ]  <-±  -  -*)  m  (  1 

AL  x  e  y  c 


Qi  A  V 

“  •  *  V  \  '?7  *  57'  • 


where  i  ,  a  ,  and  "I*  «**•  unit  vectors  in  the  x,  y,  and  s  dii-*ct  ion;. , 
k  y  z 

respectively ; 


r  =  [x  +  y  +  (x  -  h)  ] L/  ‘  , 


n 


r9  =  [x“  +  +  (x  +  h)2]1'7  , 


Ql  =  Q(t  -  r ,/v)  , 


Q2  =  Q(t  -  r2/v>  . 


The  above  Hertz  vector  has  been  converted  to  spherical  coordinates 

and  the  field  calculated  from  the  equations  for  E  and  B.  These  results 

will  be  given  in  terms  of  their  spherical  components  E^,  Fv ,  E  ,  B  ,  H,,, 

and  B.  where  E  =  E  a  +  Eran  +  E.aA  and  B  =  B  a  +  B-TaT-  +  B  ~a  . 

9  rr  00  9  rr  00  f  ^ 

ALi  9 

Lr  =  4TTc —  tsin  6(4  COS  ♦  +  n  sin  *)C-Qi»2(?*  3’  ~'2'  ~3j  *  urw  Qpo(j;  ~  )  1 


*  n  cos  n[-Q1>2(2,  3;  2,  3)  ♦  pew  Q^U;  1)1 


♦  0^(3,  4,  5)(r  -  h  cos  6)[r  n  cos  9  -  hn  +  r  sin  f>U  cos  $  +  m  sin  t )  J 


+  09(3,  4,  5)(r  +  h  cos  0)[r  n  cos  0  +  hn  -  r  sin  0(1  cos  $  ♦  m  sin  t ) J ) . 


AL 


n  s  4^7-  (cos  <(l  cos  i  *  m  sin  ♦)C-Q1  .^<2 .  3;  -3)  ♦  pro  -i'i 


♦  n  sin  0[0r2(2,  3;  ?,  3)  -  utu>?  Qx 1 ;  1)] 


-  h  sin  0[hn  -  r  sin  9(1  cos  ♦  ♦  m  sin  ♦)JlQJ(3,  4,  $)  ♦  Q..(  3,  4,  M] 


t  hnr  sin  0  cos  efQ^J,  4,  5)  -Q?(3,  4,  •. )))  . 
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4ilc 


{(l  sin  4  -  m  cos  4>)[Q1,?(2,  3;  -3)  -  pew2 


H  =  — —  (  jujye)  h  sin  0(H  sin  $ 
r  4  lit 

o 


-  m  cos  <t>)  Q1,2(2,  3;  2,  3) 


™~  ( jupr. )(  •£  sin  ♦  m  cos  4>)[r  2< 2 ,  3;  -2,  -3)  -h  cos  0  01>?(r,  3;  2,3)1. 


U  *  unY"  ,tQ1>2(2,  3;  2,  3)[rn  sin  6  +  h(i  cos  #  +  m  sin  4>)J 

o 


-  Q1>2(2,  3;  -2,  -3)  r  cos  0(1  cos  $  +  m  sin  ♦  >} 


Where 


A  Q1  % 

0  ( 1  •  +  1 )  -  —  ♦  —  , 
1  1  rl  r2 


A  1  Qs  Q{  Q: 

Qt(3,  **»  5)  -  -  (£)2  ♦  j3(H)  —  ♦  3  ,  i  = 


A  Ql  Q,  Q0  Q-> 

Q1>2(2,  3;  i  2,  ±  3)  =  j(^)  -i  +  -i  t  j(H)  t  . 

rl  rl  r2  r2 


In  addition  to  these  calculations,  the  following  has  been  per¬ 
formed  as  a  check  on  these  results.  A  Hertz  vector  was  defined  for  the 
original  dipole  in  terms  of  a  spherical  coordinate  system  which  had  the 
dipole  at  its  origin,  and  tha  field  conditions  were  calculated.  The 
same  was  performed  for  the  image  dipole  and  its  coordinate  system.  The 
two  fields  were  superimposed  and  converted  to  the  coordinate  system 
used  in  tha  first  calculations,  and  the  results  were  found  to  agree 
with  those  givon  above. 

In  the  event  that  it  is  desired  to  have  the  field  conditions  ex¬ 
pressed  in  term*,  of  a  spherical  coordinate  system  translated  in  the 
xy-piane,  the  following  result  is  needed.  Tor  example,  see  figuie  io 
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in  the  first  section.  Here  the  (x,  y,  z)  coordinates  are  translated 

by  the  vector  (x  ,  y  ,  o)  with  respect  to  the  x  ,  y  ,  z  )  coordinate 
“0  o  o  o 

system.  That  is  (xq,  yQ,  zq)  =  (x  t  xe>  y  +  yfi,  z).  The  transfor¬ 
mation  of  a  vector  A  =  a  A  +  A.  +  A,  into  a  vector  A  =  a  A 

r  r  9  9  4  4  rr 

o 

♦**0  Ao  +  %  %  is  8ivon  bv 

o  o  o  o 


A„  "  A  {(sin  0  cos  9)  sin  0  cos  4  + 

r  r  o  o 


+  A .{(cos  t)  cos  4)  sin  o  cos  a  t 

0  o  c 


+  A. {-(sin  4}  sin  0  cos  4  ♦  (cos 

4  o  o 


A  =  A  {(sin  9  cos  4)  cos  6  cos  4  + 

0  ^  r  T  o  o 


+  A.Kcos  0  cos  4)  cos  9  cos  4 
0  o  o 


+  A  {(-sin  4)  cos  0  cos  4  *  (cos 

f  o  o 


and 


(sin 

G 

sin 

4) 

sin 

0 

sin 

4 

+ 

( cor,  0 ) 

C  o 

0 

o 

(cos 

0 

sin 

4) 

sin 

e 

sin 

4 

(sin  (i) 
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o 

o 

4)  sin 

0  ! 
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4  ) 

f 

0 

o 

(sin 

0 

s 5  »-• 

4) 

COS 

0 

sin 

4 

(cos  ♦  ) 

s  i: 

o 

0 

(cos 

e 

sin 

4) 

cos 

0 

o 

sin 

*o 

* 

(sin  *•) 

■;  1  r 

4)  cos  9  sin  4  }, 

o  o 


Aa  =  A  {  -  (sin  0  cos  4)  sin  4  t  (sin  P  sin  4)  cos  4  } 
$  r  o  o 

o 

♦  A  {  -  (cos  0  cos  4)  sin  *  ♦  (cos  P  sin  4)  cos  4  } 

*'  o  o 

♦  A  ((sin  4)  sin  4  ♦  (cos  4)  cos  4  ), 

4  0  o 

where 

ro 

COS  On—  cos  0  , 

r  o 

sin  0  a  (1  -  <;£>2  ce»?  d  )i/?  . 

r  o 

sin  4  8  <r  sin  e  sin  %  -  y.)  /  (r?  -  r?  cos?  4  )i/J  , 

o  o  o  '•  o  o 

and 

c<w»  4  «  (r  sin  9  cos  4  -  X  )  /  (r?  -  r?  <os?  0  )i/?  . 

o  o  o  •  o  o 
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As  in  the  electric  dipole  calculations,  the  field  due  to  a 
m, ic  dipole  above  the  earth's  surface  may  be  calculated  by  con- 
during  the  actual  dipole  and  its  imago  radiating  in  free  space. 
Thu:.,  according  to  the  principle  of  complementary  solutions,  the 
■  i<  Id  expression;;  f- ?’  the  magnetic  dipoles  may  be  obtained  from 
f  ;.«>  .*?  of  the  ei»*f  trie  dipole  by  the  relations 

E'  =  -  ✓  1/iiC 


and 

B'  =  ^  uc  E  . 

There  is  only  one  additional  point  to  note.  Tor  an  electric 

dipole  with  orientation  (lfi,  mp,  ne>,  the  image  will  have  orientation 

(-1  ,  -m  ,  n  ),  but  for  a  magnetic  dipole  with  orientation  (1  ,  m  ,  n  ), 
ce  e  v  m  m  n 

it:;  image,  due  to  an  assumed  infinite  conducting  plane,  will  have 

orientation  (1  ,  m  ,  -n  ).  Thus,  i  >  the  above  expressions  replace 
mm  m 

1  by  -1,  n  by  -m,  and  n  bv  -n  when  these  are  associated  with  the  r^  terms. 
Model  Improvement 

The  basic  modet  related  to  the  antenna  still  contains  all  of  the 

unknowns  it  originally  h*  1  at  formulation: 

i.  The  time  distribution  and  its  associated  parameters. 

?.  The  space  distribution  and  its  parameters. 

3.  The  orientation  parameters, 

•».  The  Touricr  coefficients  of  the  waveform. 

i.  The  number  of  elements  and  their  positions  required  to 
build  up  the  "improved  model". 

It  is  suspected  that  the  filth  unknown  will  be  strongly  affected  bv 
the  first  four,  or  possibly  absorbed  into  them. 

The  work  to  be  done  in  model  improvement  consists  of  using  the 
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antenna  <iata  to  learn  the  unknowns  1 is  tea  above.  These  unknowns 
contain  the  parameters  of  parametric  pattern  recognition.  At  this 
point  only  the  first  of  the  five  has  been  considered  in  any  detail. 
The  remainder  of  this  section  will  be  devoted  to  results  obtained 
that  concern  this  time  distribution. 

The  data  used  to  investigate  the  time  distribution  was  obtained, 
from  an  isolated  cell  in  New  Mexico.  Only  the  sferic  bursts  which 
correspond  to  visually  confirmed  lightning  strokes  have  been  utilise.! 
at  the  present  time.  Except  for  a  very  few  strokes  which  .seemed  to 
trigger  several  succeeding  strokes,  the  data  appeared  to  satisfy  the 
basic  postulates  of  a  Poisson  process.  Thus  the  Poisson  process  was 
chosen  as  a  likely  candidate. 

The  distribution  function,  F^.(t),  for  the  time  between  strop 

may  be  obtained  as  follows: 

P(T  >  t)  =  Prob.  [no  strokes  in  the  interval  (o,  T)] 

-At 

=  e 

FT(t)  =  F(T  <  t)  =  1  -  P(T  >_  t)  =  l  -  e'Xt  . 

The  above  equation  may  be  rearranged  to  yield 
-At  =  ln[l  -  FT(t)3  . 

A  magnitude  scaled  sample  distribution,  F^(t),  has  been  obtained 
from  the  data.  The  scale  factor  is  214  since  there  were  214  obser¬ 
vations.  Thus  if 

.  F®(t)  =  214(1  -  e'At)  , 

-At  =  -  In  214  +  ln( 214  -  r®(t)  )  , 


then 


2 

and  a  plot  of  (214  -  F\p(t)  verses  t  on  semi-logarithmic  paper  should 
yield  a  straight  line  with  slope  -X.  Those  results  are  shown  for 
*  -  .oiiQV  in  figure  12.  figure  13  shows  l'^(t)  verses  t  in  dots  and 
lli<*  (unction  214(1  -  e  as  a  continuous  lino,  from  these  graphs, 

it  appears  that  the  function  r^.(t)  =  1  -  e  ‘  is  a  reasonable 
estimate  for  the  distribution  function  of  this  particular  storm. 

Consider  now  the  corresponding  density  function 


fT(t) 


dr^(t) 

t 


for  the  assumed  value  X  =  .0682,  the  scaled  density  function  in 


f^(t)  =  214  Xe"U  =  14.6  e~-0682t  . 

Tli ir.  function  is  compared  with  the  empirical  density  function  f^O) 
in  figure  14  of  this  section.  Consideration  of  the  empirical  density 
function  suggests  that  it  may  he  other  than  Poisson  with  uniform  A. 

One  form  suggested  by  the  shape  of  the  empirical  function  is  a 
Rayleigh  distribution.  However,  this  is  in  total  disagreement  with 
the  sample  distribution.  Another  possibility  is  that  the  rate  is 
Poisson  distributed,  but  X  is  not  constant. 

This  summarizes  the  work  done  to  present.  The  following  section 
will  mention  briefly  the  intended  method  of  classification  of  the  storm 
and  the  "criterion  of  goodness"  for  the  model. 

Pattern  Recognition 

How  the  model  is  to  be  used  in  classification  can  not  be  exactly 
specified  until  further  work  is  done.  It  may  be  possible  to  classify 
a  storm  by  considering  only  changes  in  the  model's  parameters,  or  the  mo. 
may  require  a  change  in  form  as  the  storm  progresses.  In  any  event,  it 
is  claimed  for  the  purposes  of  this  phase  that  any  model  which  yields 


figure  14  C  (d  versus  time 


o  "noticeable"  difference  a.;  the  storm  chanp.es  will  be  a  "good"  model. 

R.  NONPARAMETRIC  PATTERN  RECOGNITION 
hit  reduction 

Two  nonparametric  methods  of  pattern  recognition,  threshold 
logic,  and  nearest  neighbor  rules  have  been  investigated  and  programmed 
on  t ho  computer.  These  methods  are  described  briefly  in  following 
•  lions.  Both  of  these  methods  require  changing  the  form  of  the 
recorded  data  to  the  form  required  by  the  computer  (data  reduction) 
and  selecting  from  the  many  possible  measurements  those  that  are  to 
be  used  (measurement  selection)  in  pattern  recognition. 

Although  the  ultimate  goal  of  pattern  recognition  is  to  identify 
meteorological  phenomenon,  an  intermediate  project  is  to  classify 
lightning  discharges  into  two  categories:  Cloud  to  cloud  and  cloud  to 
ground.  This  project  was  selected  in  order  to  test  the  pattern  rec¬ 
ognition  programs  with  presently  available  data. 

Threshold  Logic  Unit 

The  threshold  logic  unit  (TLC)  is  a  linear  trainable  pattern 

classifying  scheme  and  can  be  used  to  classify  patterns  into  one  of 

two  categories  (Example:  Haii-no  Hail,  Tornado-no  Tornado).  The  TLU 

work:;  well  if  the  patterns  belonging  to  the  two  categories  are  linearly 

separable.  A  set  of  patterns  in  the  n  dimensional  space  <E  )  are 

linearly  separable  if  there  exists  a  hyperplane  in  such  that  all 

the  patterns  b<  onging  to  one  category  lie  on  one  side  of  the  hyper- 

plane.  Tbs  equation  of  the  hyperplane  can  be  written  as.  g(X)  = 

w.X.  ♦  w-X  ♦...  ♦  u>  X  ♦  w  ,  *  0,  where  the  "weights"  w,  ,  u>.,...w  ,  w 
1  i  *  *  n  n  nti  l  i  n  nt  t 

ai<*  such  that 

p.(X)  *■  0  If  X  belongs  to  category  1  and 

— *  — » 

g(X)  <0  If  X  belong*:  to  category  II. 


j'.(X)  is  called  the  "discriminating  function"  and  can  be  easily 
implemented  according  to  figure  15. 


FIGURE  15.  THRESHOLD  LOGIC  UNIT 

The  response  of  the  threshold  element  is  +  1  if  g()0  >  0,  and  is  -  l 
if  g(X)  <  0.  The  pattern  is  classified  as  category  I  or  II  depending 
on  whether  the  response  is  +  1  or  -  1. 

If  one  has  a  complete  knowledge  about  all  the  patterns  to  be 
classified,  then  the  weights  ,w2  ’ "  ’  *  ,UJn+l  can  determined  analyt  i¬ 
cally.  However,  in  most  of  the  practical  cases,  only  a  few  sample 
patterns  are  available.  In  such  cases,  these  weights  are  determined 
by  trial  and  error  based  on  the  sample  patterns.  The  iterative  method 
of  Arriving  at  the  final  values  for  the  weights  is  called  "training 
the  TLU" ,  and  the  patterns  used  in  the  process  are  referred  to  as  the 
"training  set."  Training  the  TLU  is  done  on  a  computer,  and  a  computer 


program  has  been  written  and  checked  out.  For  further  description  of 
TLU's  see  reference  2. 

Nearest  Neighbor  Rule 

The  pattern  classifying  scheme  based  on  the  nearest  neighbor 
rule  is  wcil  suited  for  use  in  situations  where  the  patterns  belonging 
t'<  i  particular  category  tend  to  cluster.  (Figure  16). 
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Figure  lb.  NEAREST  NEIGHBOR  RULE 


As  the  name  implies,  the  machine  classifies  the  pattern  X  as  coming 
from  category  j,  if  th;  training  pattern  nearest  to  X  is  of  categorv 
).  A  measure  of  nearness  is  the  Eueledian  distance  given  by 


d(x"",  pj)  =  *'  (x  -  p!  )  •  oT  -  f  j  ) 


k  |  ^ 

win  re  d ( X ,  f.)  r  the  distance  between  X  and  the  i  training  rut  tens 
t  h  . 

in  the  1  category.  Tito  distance  between  X  .aid  the  nearest  pittern 

.  .  t  h  ,  , 

of  i  category  is  given  by 


d(j)  *  min  C  <1<  X ,  P.1))  , 
i  1 


and  this  ir.  celled  D(i>,  t  tie  distance  between  X  a»  d  tin'  j‘  '  -.it  ecu  ,  . 


X  is  classified  as  category  j  if 


L)(  j)  <  D(k) ,  k  =  1,  2  .  .  .  P.  ;  k  i  j 


where  R  is  the  total  number  of  categories. 

For  further  description  of  nearest  neighbor  rules  see  reference  3. 
A  computer  program  tor  nearest  neighbor  pattern  classification  has  boon 
written  and  cheeked  out. 

Data  Reduction 
Narrow  Band  Recording 

The  present  set-up  for  digitising  the  analog  data  from  the 
narrow  band  recording  was  outlined  in  the  quarterly  repoi’t  dated  June 
30,  1068.  On  this  set-up,  the  analog  data  is  digitised  using  an 
F.IXIO-761  A-D  converter,  and  the  digital  sampler,  are  fed  into  the  core 
memory  of  an  IBM  1020  computer.  After  the  core  is  filled,  its  contents 
are  transferred  to  an  IBM  1311  disk  unit  and  stored  there  for  further 
processing. 

The  analog  to  digital  conversion  is  done  on  one  track  at  a  time, 
at  the  rate  of  600  samples  per  second.  In  order  to  be  able  to  look 
at  identical  sections  of  data  on  all  the  12  tracks  of  data,  a  common 
time  reference  must  be  established.  This  time  reference  is  estab¬ 
lished  by  masking  the  tape  for  a  few  inches.  The  masking  produces 
•/.oro  analog  output  voltage  on  all  the  twelve  channels,  and  this  voltage 
appears  as  a  string  of  zeros  in  the  digital  samples.  This  string  of 
zeros  can  be  easily  identified  and  used  as  a  time  reference  for  locatin, 
sections  of  digital  data  on  different  channels  corresponding  to  a 
particular  interval  of  time. 

Computer  programs  have  been  developed  for  further  processing  of 
the  digital  data  from  the  narrow  band  recording.  Programs  arc  now  i*. 


use  for  calculating  the  average  value  of  the  signal,  the  amplitude 
distribution,  and  the  sferic  count,  and  the  correlation  between  the 
signals  in  the  horizontal  and  vertical  channels  for  each  frequency. 

At  present,  the  digital  data  stored  on  IBM  1311  disks  have  to  he 
transferred  to  punched  cards  before  further  processing  can  be  done. 
This  prevents  a  detailed  analysis  of  large  segments  of  narrow  band 
data.  However,  new  computing  facilities  are  expected  to  be  available 
at  the  end  of  fall  1968  (see  section  on  data  collection).  This 
includes  an  IBM  360  computer  which  will  accept  input  data  directly 
from  dick  storage,  thus  eliminating  the  need  for  transferring  the 
data  to  punched  cards. 

Wide  Band  Recording 

The  sampling  rate  required  for  analog  to  digital  (A-D)  con¬ 
version  of  the  wide  band  recording  is  about  600,000  samples  per 
second.  The  present  A-D  conversion  equipment  can  handle  a  maximum 
rate  of  500  samples  per  second.  The  wide  band  recording  cannot  be 
digitised  at  this  rate  without  loosing  much  of  the  high  frequency 
components  of  the  signal.  However,  by  expanding  the  time  base  of  the 
wide  band  recording  by  a  ratio  of  1:1024,  the  effective  sampling  rate 
can  be  increased  to  512,000  samples  per  second.  At  this  sampling 
rate,  it  is  possible  to  recover  the  components  of  the  signal  up  to  a 
frequency  of  250  KHz. 

A  time  base  expansion  of  1:1024  is  obtained  in  two  steps: 

1.  The  original  tape  containing  the  recording  taken  at  60  ips 
is  played  back  at  a  speed  of  1  7/8  ips  and  recorded  on  a 
new  tape  at  60  ips.  The  new  tape  now  has  a  time  base 
expansion  of  1:32. 

2.  The  new  tape  when  replayed  at  a  speed  of  1  7/8  ips 
expands  the  time  base  by  another  factor  of  1:32,  giving 
a  total  time  base  expansion  of  1:1024  (322). 

The  digitising  of  wide  band  recording  is  in  the  final  testing 


stages.  Programs  have  already  been  developed  to  Fast  Fourier  transform 
the  digital  data.  Fast  Fourier  transforming  is  in  effect  an  averaging; 
hence,  it  reduces  the  number  of  data  points  in  the  wide  band  data. 

Also,  this  gives  a  way  of  comparing  the  power  spectrum  so  obtained 
with  the  narrow  band  mean  square  averaging. 

The  wide  band  recorder  has  a  frequency  response  (±  3db)  of 
300-300,000  Hz  at  a  recording  speed  of  60  ips.  Hence,  the  signal 
coming  off  the  wide  band  recording  will  have  high  frequency  components 
up  to  and  beyond  300  KHz.  Because  of  the  slower  sampling  rate  of 
51?, 000  per  second,  the  power  density  spectrum  of  the  signal  obtained 
from  these  digital  samples  will  have  some  error  due  to  "aliasing." 
However  an  inspection  of  the  narrow  band  (mean  square  averaged) 
recording  of  the  signal  indicates  that  the  energy  content  is  small  at 
high  frequencies  beyond  250  KHz.  Therefore,  error  due  t:o  aliasing 
can  be  neglected. 

Measurement  Selection 

The  digital  samples  obtained  from  the  narrow  band  and  wide 
band  recordings  constitute  the  data  for  the  pattern  recognition  scheme. 
The  narrow  band  comb  filter  outputs  have  been  recorded  continuously, 
while  the  wide  band  signal  has  been  recorded  for  one  minute  duration 
once  in  every  fifteen  minutes.  In  order  to  be  able  to  use  the  data 
from  both  the  recordings  to  form  pattern  vectors,  it  is  necessary  to 
take  segments  of  narrow  band  data  for  fifteen  minutes  duration  and 
combine  it  with  the  corresponding  wide  band  data  for  one  minute  duration. 
The  combined  data  for  this  fifteen  minute  period  is  used  to  form  one 
pattern  vector. 

At  the  ideal  sampling  rate,  there  are  36  million  digital  samples 
per  minute  from  the  wide  band  recording  and  about  360,000  samples  per 
minute  from  the  narrow  band  recording.  It  is  not  possible  to  use  such 
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a  large  number  of  samples  to  describe  patterns.  Instead,  a  reasonable 
number  of  measurements  (variables )  characterising  the  data  will  be  used 
to  form  the  pattern  vectors. 

A  set  of  44  measurements  have  been  chosen  tentatively  to  de¬ 
scribe  the  characteristics  of  the  data  taken  over  a  period  of  fifteen 

minutes.  Each  of  these  variables  is  quantized  into  five  levels,  re- 
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suiting  in  a  maximum  of  5  distinct  patterns.  Quantizing  is  done  to 
save  memory  storage  on  the  computer  and  to  increase  the  speed  of 
computation.  The  nature  of  these  measurements  and  their  levels  of 
quantizing  are  explained  below. 

Measurements  From  Wide  Larid  Recording 

Because  of  the  large  number  of  data  samples,  preprocessing  is 
necessary  on  the  wide  band  data  before  any  measurements  can  be  at¬ 
tempted.  Fast  Fourier  transforming  will  be  a  preprocessing  on  wide 
ba  d  data,  and  a  set  of  14  measurements  have  been  chosen  to  char¬ 
acterize  the  power  density  spectrum  obtained  by  Fast  Fourier  trans¬ 
forming.  These  measurements  describe  the  change  in  power  density 
between  two  successive  bursts  of  wide  band  recording  at  the  following 
frequencies:  DC  (recovered  from  wide  band  T.M.  recording),  10  KHz, 

25  KHz,  50  KHz,  .  .  .  300  KHz.  The  differences  in  power  density, 
instead  of  the  actual  values,  are  taken  because  of  two  advantages: 

1.  By  taking  the  difference  in  power  density  corresponding 

to  two  successive  segments  of  data,  some  information  about 
the  past  history  of  the  signal  is  maintained. 

2.  Any  bias  or  error  in  the  measurements  tend  to  cancel  out. 

Each  of  these  14  measurements  are  quantized  into  five  levels  as 


indicated  below. 

Amount  of  change  in  power  density  Quantized  value 

Above  50%  decrease  compared  to  previous  value  .  1 
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50-5%  decrease  compared  to  previous  value  .  2 

With  in  5%  of  the  previous  value . 3 

5-50%  increase  over  previous  value  .  4 

Above  50%  increase  over  previous  value  .  .  .  .  , . .  5 

rteasuretnents  From  Narrow  Band  Recording 


A  set  of  30  measurements  are  taken  to  represent  the  characteristics 
of  the  data  coming  from  the  12  channels  of  the  narrow  band  recording. 
Those  measurements  contain  information  about  the  sferic  counts  on 
each  channel,  the  amplitude  distribution,  and  the  correlation  between 
the  signals  in  the  horizontal  and  vertical  channels. 

Out  os  these  33  measurements,  a  set  of  12  are  chosen  to  measure 
the  sferic  counts  on  the  12  channels  of  data.  The  sferic  count  is 
taken  with  respsct  to  a  threshold  level,  which  is  set  slightly  higher 
than  the  average  value  of  the  signal  in  that  channel. 

<VP  -  VA) 

VT  =  t  ' - ft .  -  Average  value 

Vp  -  Peak  value 
Vp  -  Threshold  value 

The  sferic  counts  will  be  different  on  each  channel,  and  the  values 
of  the  sferic  counts  are  quantised  with  respect  to  the  average  count 
over  the  12  channels. 


Actual  Values  of  the  Sferic  Count  Quantized  value 

Less  than  50%  of  the  .verage  value  . . 1 

Between  50-95%  of  the  average  value  . 

Between  95-105%  of  the  average  value  ....  .  .3 

Between  105-150%  of  the  average  value  ........  .  u 

Over  150%  of  the  average  value  . . 5 


Characteristics  of  the  amplitude  distribution  of  the  data  are 
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presented  by  a  set  of  12  measurements.  Each  measurement  represents 
the  percentage  of  total  time  during  which  the  signal  in  a  particular 
channel  is  above  the  threshold  value  set  for  that  channel.  The  actual 
values  of  the  measurements  on  the  amplitude  distribution  are  quantized 
as  follows. 


Actual  value  of  time 

he low  20%  ...  . 

Between  20-40%  . 

Between  40-60%  ..... 

Between  60-80%  . 

Above  80%  . 


Quantized  value 

. 1 

. 2 

. 3 

. 4 

. 5 


Six  measurements  are  taken  to  represent  the  correlation  between, 
signals  in  the  horizontal  and  vertical  channels  at  frequencies  10  KHz, 
50  KHz,  100  KHz,  150  KHz,  200  KHz  and  250  KHz.  The  actual  values  of 
correlation  u :*e  quantized  as  follows: 

Actual  Valuer  of  Correlation  Quant iz'ed  value 


Below  0.2 . 1 

Between  0.2  to  0.4  . . . . . . 2 

Between  0.4  to  0.6  o 

Between  0  6  to  0.8  .........  ...............  4 

Above  0.8.... . 5 


At  this  stage ,  it  is  not  known  if  these  measurement  will  be 
adequate  to  transform  the  data  into  recognizable  patterns.  However, 
•liter  the  initial  training  of  the  p.-ittem  recognition  scheme,  new 
measurements  will  he  added  if  nec.^oary.  Also,  some  meas'-ementr. 
may  be  omitted  If  the  weights  on  the  TLU  corresponding  to  these  meant;- 


5b 


merits  are  very  small  compared  to  others.  Tor  further  discussion  on 
measurement  selection  see  reference  4. 

An  Intermediate  Pattern  Classification  Test 

As  mentioned  earlier  in  this  section,  two  nonparametric  methods 
of  pattern  recognition— threshold  logic  unit  and  nearest  neighbor 
rule-will  be  used  to  classify  severe  weatl  er  patterns.  Two  major 
steps  involved  in  implementing  these  pattern  recognition  shcemes  are 

1.  Taking  measurements  from  the  data  to  form  the  training 
patterns. 

2.  Training  the  computer  to  learn  to  classify  the  patterns  in 
the  training  set. 

Computer  programs  have  been  developed  for  obtaining  measurements  from 
the  data  to  form  pattern  vectors  and  to  classify  them  using  the  two 
nonparametric  methods.  Attempts  are  now  being  rade  to  classify  lightning 
discharges  into  two  categories  (cloud  to  cloud  or  cloud  to  ground) 
vr.ing  these  programs.  A  brief  outline  of  the  procedure  is  given  below. 

The  edge  tracks  of  the  narrow  band  recordings  taken  during  some 
recent  storms  contain  K  ief  descriptions  about  the  nature  of  the 
lightning  discharges  based  on  visual  observation  from  the  aircraft. 

This  description  includes  such  details  as  the  nature  of  the  discharge 
and  if  the  discharge  occurred  close  to  the  aircraft.  Py  correlating 
with  this  information,  any  particular  burst  of  activity  in  the  narrow 
band  channels  can  be  identified  as  due  to  a  cloud  to  cloud  discharge 
or  cloud  to  ground  discharge.  Sections  of  narrow  band  data  corre¬ 
sponding  to  individual  discharges  have  been  digitised  and  plotted 
using  an  IBM  1620  computer  and  Cal<  ump  plotter. 

Figure  17  and  18  show  the  digital  plots  of  the  narrow  band  signal 
due  to  two  cloud  to  cloud  discharges.  The  signals  in  the  horizontal 
and  vertical  channels  at  three  frequencies  (10  KHz,  SO  KHz,  and  ISO  Khz) 
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TOP  -SIGNAL  IN  10  Khz  HORIZONTAL  CHANNEL 
BOTTOM  -  SIGNAL  IN  10  Khz  VERTICAL  CHANNEL 

TIME  SCALE:  l"  “  100  MILLISECONDS 


FIGURE  17a.  CLOUD  TO  CLOUD  DISCHARGE  I. 


i 
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TOP  -  SIGNAL  IN 
BOTTOM-  SIGNAL  IN 


50  Khz  HORIZONTAL  CHANNEL 
50  Khz  VERTICAL  CHANNEL 


FIGURE  17b.  CLOUD  70  CLOUD  DISCHARGE  I. 
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TOP  -  SIGNAL  IN  150  Khz  HORIZONTAL  CHANNEL 
BOTTOM- SIGNAL  IN  150  Khz  VERTICAL  CHANNEL 

TIMESCALE:  l"  «  100  MILLISECONDS 


FIGURE  17c  CLOUD  TO  CLOUD  DISCHARGE  I 


TOP  -SIGNAL  IN  10  Khz  HORIZONTAL  CHANNEL 
BOTTOM- SIGNAL  IN  10  Khz  VERTICAL  CHANNEL 


FIGURE  I80  CLOUD  TO  CLOUD  DISCHARGE  II 
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TOP  -SIGNAL  IN  50  Khz  HORIZONTAL  CHANNEL 
BOTTOM -SIGNAL  IN  50  Khz  VERTICAL  CHANNEL 

TIME  SCALE  l“  ■  100  MILLISECONDS 


FIGURE  ClOUO  TO  CLOUD  DISCHARGE  n 
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TOP  -  SIGNAL  IN  150  Kh2,  HORIZONTAL  CHANNEL 
BOTTOM-  SIGNAL  IN  150  Knz,  VERTICAL  CHANNEL 

TIME  SCALE:  f  *  100  MILLISECONDS 


FIGURE  UK  CLOUD  TO  CLOUD  DISCHARGE  II. 


e>2 

are  shewn  in  these  figures.  Figure  19  shows  the  plot  of  a  cloud  to 
ground  stroke.  At  each  frequency  the  gain  is  the  same  for  the  hori¬ 
zontal  and  vertical  channels.  An  inspection  of  these  plots  leads  tc 
the  belief  that  there  is  nn.i'e  activity  in  the  vertical  channels  during 
a  cloud  to  ground  discharge  compared  to  a  cloud  to  cloud  discharge. 

(It  is  hoped  that  the  computer  win  recognize  this  difference,  through 
our  measurements) . 

A  set  of  30  measurements  taken  from  the  narrow  band  data  is  used 
to  fozmi  training  patterns.  This  set  of  measurements  is  the  same  as 
the  one  outlined  in  section  1I-D  (on  measurement  selection).  Training 
patterns  for  each  category  are  now  being  formed.  Pattern  recognition 
will  be  attempted  after  a  reasonable  number  of  training  patterns  nave 
been  obtained.  Then  the  trained  computer  will  be  tested  on  other 
samples . 

The  data,  the  measurements,  and  the  techniques  used  here  +o 
classify  lightning  discharges  are  similar  to  the  ones  that  will  be 
used  later  to  classify  severe  weather  patterns.  Therefore,  it  is 
hoped  that  the  experience  gained  through  this  intermediate  attempt 
will  be  helpful  when  classifying  meteorological  phenomenon. 
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TOP  -  SIGNAL  IN  10  Khz 

HORIZONTAL  CHANNEL 

BOTTOM  -  SIGNAL  IN  10  Khz 

VERTICAL  CHANNEL 

TIME  SCALE:  l"*  100  MILLISECOND 


FIGURE  19a  CLOUD  TO  GROUND  STROKE 
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TOP  -  SIGNAL  IN  50  Khz 

HORIZONTAL  CHANNEL 

BOTTOM  -  SIGNAL  IN  50  Khz 

VERTICAL  CHANNEL 

TIME  SCALE:  lH«  100  MILLISECONDS 


TO  GROUND  STROKE. 
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TOP  -  SIGNAL  IN  150  Khz 

HORIZONTAL  CHANNEL 

BOTTOM  -  SIGNAL  IN  150  Khz 
VERTICAL  CHANNEL 

TIME  SCALE:  lN  *  100  MILLISECONDS 


FIGURE  I9c  CLOUD  TO  GROUND  STROKE. 
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V.  DrtTA  ANALYSIS 

A.  USE  OF  EXISTING  FACILITIES 

An  IBM  7040  digital  computer  is  presently  available  for  processing 
data  recorded  in  the  aircraft.  To  date,  the  7040  har.  been  used  pri¬ 
marily  for  program  development  and  checkout;  also,  some  preliminary  runs 
have  been  made  on  actual  atmospheric  data. 

Within  the  month,  the  7040  machine  will  be  replaced  by  an  IBM 
360/50  computer. 

Before  either  of  the  computers  can  be  used,  the  raw  data,  recorded 
on  magnetic  tape  in  analog  form,  must  be  digitized  and  put  on  a  media 
acceptable  as  input  to  the  digital  computer.  In  addition,  the  higher 
frequency  data  must  be  "time-scaled"  by  slowing  down  the  tape  machine 
speed  on  playback. 

Selected  portions  of  low-frequency  data  are  presently  being 
digitized  and  the  results  punched  into  cards  through  the  use  of  an 
existing  slew-speed  system.  This  system  is  extremely  inadequate  as 
far  as  the  preparation  and  storage  of  long  runs  and  many  runs  of 
data  are  concerned;  however,  it  is  serving  to  furnish  realistic  data 
for  the  testing  of  the  various  programs. 

Briefly,  the  characteristics  and  details  of  operation  of  the 
sy.  tem  arc  shown  in  Figure  20,  The  data  tape  is  played  back  at  the 
ground  laboratory,  and  the  rjgnal  is  conditioned  in  a  preamp  to  a 
^standard  range  of  o-*»  volt:-..  The  signal  is  then  fed  into  an  analog 
to  digital  convorU-r  which  converts  the  samples  to  BCD  digital  form. 

Each  sample  consists  of  throe  decimal  digits. 


DIGITS 

MAX  Sample  rate  =  500  Sar..;;les/sec 

MAX  Duration  of  one  samp*"  run  =  6  sec 

TIME  Required  to  dump  core  to  disk  =  1/20  sec 

Figure  20.  PRESENT  DIGITAL  DATA  PREPARATION  EQUIPMENT 

The  sampling  process  continues  until  the  core  memory  of  the  1620 
computer,  shown  in  Figure  20,  is  effectively  filled.  The  computer 
then  dumps  its  core  memory  and  then  proceeds  to  resume  the  sampling 
process. 

The  speed  and  volume  handling  capabilities  of  the  system  are 
primarily  limited  by  the  speed  and  core  memory  size  of  the  1620  com¬ 
puter.  The  system  can  convert  data  at  a  maximum  rate  of  500  samples 
per  second  for  an  uninterrupted  duration  of  3000  samples.  At  this 
time,  the  core  memory  must  be  dumped  v.ito  the  disk,  and  this  requires 
approximately  0.05  second. 

Later,  the  data  on  the  disk  is  punched  into  cards,  20  samples 
per  card.  Obviously,  the  digitising  of  even  a  moderate  amount  of 
data  results  in  an  exorbitant  numl-r  of  cards. 

Now  Equipment 

In  ordar  to  seriously  pursue  the  preparation  of  the  data  obtained 
from  the  flight  recorder-;,  new  equl:.-sent  must  be  secured.  In  brief, 
thin  in  because  sampling  rates  must  go  an  high  as  20,000  samples  per 
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second  for  extended  runs  involving  up  to  approximately  2  million 
samples.  So,  not  only  must  this  supporting  equipment  possess  a  much 
higher  speed  capability  than  the  existing  equipment,  but  it  also 
must  have  a  much  larger  and  higher  density  storage  media.  This 
dictates  digital  magnetic  tape  as  the  output  storage  media  of  the 
system;  and,  the  form  of  the  digitized  data  on  this  tape  must  be 
compatible  with  the  IBM  ju0/50  computing  system. 

Furthermore,  since  the  360/50  system  will  not  have  a  graphical 
plotting  feature,  it  is  virtually  imperative  that  any  proposed  new 
equipment  should  include  this  capability.  A  Caicomp  plotter  is 
presently  available,  so  the  equipment  will  only  need  the  necessary 
interfacing  and  control. 

Some  of  the  recorded  flight  data  represents  analog  atmospheric 
energy  with  information  frequencies  as  high  as  250  kilohertz,  approx¬ 
imately.  The  tape  with  this  data  can  be  sloweu  down  by  a  ratio  of 
32  to  1  upon  playback,  thus  effectively  frequency  scaling  the  data 
by  the  same  ratio.  The  scaled  data  can  then  be  said  to  contain 
frequencies  of  up  to  250/32  =  8  KHz. 

In  order  to  adequately  represent  this  scaled  data  in  digital 
form,  a  minimum  sampling  rate  of  approximately  20,000  samples  per 
second  is  dictated. 

A  typical  burst  of  atmospheric  signal  will  require  a  total 
number  of  samples  in  the  range  of  1  to  2  million  without  any  gaps 
in  the  data  vector  in  order  to  mat>*  possible  the  desired  analysis 


procedures. 


A  number  of  digital  conversion  systems  have  been  investigated 
and  quotations  were  obtained  on  several.  The  system  that  combines 
an  acceptable  cost  with  the  desired  specifications  of  sample  speed, 
tape  input,  tape  output  (9  channel,  IBM),  plotter  interface,  and  a 
flexible  method  of  data  editing  control,  is  the  Hewlett-Packard  2115A 
system  equipped  with  an  HP3030  tape  machine  and  an  EEC0761A  analog 
to  digital  converter. 

Steps  have  already  been  ini  tinted  to  purchase  this  system;  once 
the  order  is  received  by  the  Hewlett-Packard  Company,  delivery  is 
expected  to  be  12-14  weeks. 

In  summary,  the  new  system  being  proposed  will  accept  analog 
data  and  digitize  it  at  rates  of  up  to  at  least  20,900  samples  per 
second.  The  output  will  be  put  on  a  nine-track  magnetic  tape  in  a 
form  compatable  vith  the  IBM  360/50  system.  This  sampling  task  can 
continue  until  the  output  tape  is  filled.  The  output  tape  will 
accommodate  approximately  10  million  samples  if  desired. 

Output  from  the  analysis  programs  (also  on  nine- track  magnetic 
ta;>e)  run  on  the  3C0/50  computer  can  be  read  into  the  HP2115A  system 
and  plotted  on  a  Cuicomp  plotter. 
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B.  DIGITAL  DATA  RfDUCTION  AND  TR' FORMATION 

The  exact  data  reduction  and  analysis  processes  to  be  applied 
to  the  collected  data  are,  as  yet,  not  known.  However,  it  seems 
that:  Fourier  Analysis  will  undoubtedly  be  of  primary  interest  in 
study  of  the  wide-band  signal  recordings.  With  the  wide-band  analog 
data  sampled  and  recorded  in  digital  form  (suitable  for  direct  input 
to  the  University's  IBM  360/50  computer),  the  Finite  Discrete  Fourier 
Transform  is  applicable  for  such  analysis.  The  Cooley- Tukey  algorithm 
[b]  provides  an  efficient  high-speed  implementation  for  the  necessary 
processing. 

However,  because  of  the  large  bandwidth  of  the  signal  (0-25C  KHz), 
meaningful  analysis  of  data-tracks  only  one  or  two  secon  ’ength 

will  involve  large  streams  of  data.  Samples  of  one  or  two  million 
observations  may,  on  occasion,  be  of  interest.  The  direct,  in-core 
manipulation  of  such  large  data-vectors  is  not  possible  in  the  avail¬ 
able  digital  computer.  Therefore,  development  of  a  computer  program 
( in  Fortran  IV)  has  been  carried  out  to  perform  the  Fast  Fourier 
Transform  (FFT)  on  large  data-vectors. 

A  highly-efficient  program  for  performing  the  ITT  on  in-core 
data-vectors  is  available  from  the  SHARE  library  [6],  This  program 
is  used  as  the  nuclc  is  for  the  large-vector  program  which  has  been 
developed.  As  outlined  by  Gentleman  and  Rande  [7}  and  Singelton  [8], 
the  TFT  can  Lv  p^rtormed  on  subsections  of  a  largo  data-vector  by  a 
series  of  pre- trans form  operation*  cn  the  segments.  The  number  of 
such  pre -trans from  phases  is  dependent  on  the  number  of  subsections 
iv' iuired  —  each  subsection  finally  being  submitted  to  the  FFT  program 


in-core. 


For  example,  suppose  that  t  -  larges*-  data-segment  that  can  be 
handled  in-core  with  the  FFT  prognm  is  4096  samples.  A  data- vector 
of  65,536  samples  can  be  handled  in  16  subsections.  This  would  require 

4 

four  pre-transform  phases  (s  =16)  followed  by  submission  of  each  of 
the  16  subsections  to  the  rrT  program. 

Follow. -ig  the  16  FFT  operations,  a  number  of  post-sorting  phases 
are  required  (again  4  for  the  example  16  sections)  to  re-order  the 
data  into  proper  sequence.  The  program  to  implement  this  has  been 
thoroughly  tested  on  the  existing  computer  facility,  an  IBM  7040. 

The  pre-transform  and  post-sorting  is  accomplished  using  four  scratch 
tapes  for  temporary  storage. 

The  360/50  installation  is  due  for  activation  in  early  November, 
and  will  have  available  the  IBM  2314  Disk  Storage  Facility  (8  cor- 
rdinated  11-disk  storage  units).  It  is  planned  to  begin  immediate 
prograrh  modification  to  utilize  the  disk  storage  faciltiy  for'the 

V  % 

pre-transform  and  post-sort  operations.  With  256kH  of  core  storage 
available  and  the  overlap  of  1/0  and  CPU  functions,  it  is  anticipated 
that  an  appreciable  improvement  in  total  execution  time  can  be 
obtained.  The  program  structure  places  no  upper  limit  on  the  data- 
vector  size  --  the  increase  of  execution  time  foi  very  large  vectors 
will  determine  the  allowable  sample  size. 

If  the  Fourier  Analysis  of  digitized  data  is  eventually  determined 
to  he  of  critical  importance,  and  »he  required  execution-time  is  a 
limiting  factor,  a  number  of  alternatives  are  available.  Hardware 
implementation  of  the  list  Fourier  Transform  bv  array-processor 
techniques  offer,  great  promise.  A  vector  of  1324  samples  can  be 


tr.iri3 formed  by  3uch  "black-boxes"  in  tiroes  of  the  order  of  40  milli¬ 
seconds.  The  time  increases  in  the  order  of  N  log^N,  where  N  is 
the  sample  size.  These  "black-boxes"are  currently  available  for 
direct  addition  to  the  IBM  360  machines  in  the  form  of  an  additional 
I/O  channel.  They  operate  in  the  "cycle-stealing"  mode,  and  allow 
normal  CPU  operation  once  they  are  commanded  to  operate.  Coupled 
with  graphics-display  devices,  they  offer  the  promise  of  "almost  al 
time"  Fourier  Analysis  and  spectral  display  for  quite  large  data- 


VI.  CLOUD  MODELING 


A  preliminary  study  of  the  basic  hydrodynamic  equations  of 
iwo-dimenr.ional  vortex  flow  as  they  appl>  to  an  idealized  tornado 
funnel  was  made  in  order  to  ascertain  whether  a  highly  restricted 
mathematical  mode’  would  produce  reasonably  -  ;curate  predictions. 

Tin?  results  obtained  indicate  th.it  this  preliminary  effort  was 
successful. 

The  results  of  the  tornado  investigation  are  .'hown  on  the 
following  pages  of  this  section. 

References  S  through  14,  Appendix  A,  were  used  in  the  in¬ 
vestigation  of  the  idealized  tornado  funnel. 

A.  j ASSUMPTIONS 

1.  The  tornado  energy  is  supplied  by  water  condensation  in  the 
funnel. 

'N  2.  The  rising  heated  air  in  the  funnel  is  equivalent  to  a. 
vertical-cone-sink  for  the  inward  spiralling  air. 

3.  The  radial  and  tangential  air  velocities  do  not  depend 
on  time,  and  de  end  only  weakly  on  height. 

4.  The  air  outside  the  funnel  is  approximately  incompressible. 

5.  The  heating  of  the  air  by  its  own  sheer- friction  is  negligible. 

B.  GENERAL  EQUATIONS 

1.  Navier- Stokes  Equations 

V  •  (pV)  ♦  Jfi.  *  0 

!■£.  ♦  (V*V)  V  *  ulLEll!  *  .jji  v(y.v)  ♦  jL  V^V 
t  •  time  (see) 

0  «  local  fluid  density  (ga  cm*3) 
u  1  fluid  viscosity  (dyne  sec  e*  ) 

V  *  fluid  velocity  vector  (cm  sec"1 ) 
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q  =  external  field-pot w  ial  (erg  an  ) 


<V«7)V  =  7(V»V)/7  -  VX(TXV) 


7  V  s  7(7»V)  -  VX(TXV) 

7,  7*,  7X  *  gradient,  divergence,  curl 


2.  Specialized  Equations 


q  *  0,  steady  state,  p  =  constant 
No  vertical  component  in  V 

A  A 

j  .  j  =  unit  radial  and  tangential  vectors 
r  o 


V  s  jp  p(r)  ♦  je  v(r) 


p(V»7)V  +  7p  *  p7*V 


7*V  s  0 


u(a)  *  -s  ,  v(a)  *  v 


C.  EQUATION  REDUCTIONS 


3u  .  #u  .  }v  ,  8v  Ju  3v  . 

TT  *  1  5^  3  ZT  *  TT  *  T—  a  0 


'•v  •  °*  ? S?  <r“!  *  ?lf  ■  s 


ru*  constant  •  -as 


£(¥•7)73 


au  v" 


r  Jr  r 


■  ?  h  «'  TF>  -  7 


t <»•?)»),  ■  5.JJ*  ~ 


■  ?  I?  «*■  S?»  •  T 


t 
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Limiting  case  of  zero  viscosity 

m  ■*  0  means  k  ■*  <• 

r  Tf  *  v  1  0  1?  (rv) 

This  is  satisfied  by  the  equation 
v  =  aw/r  of  a  perfect  vortex. 


Limiting  case  of  crucial  sink  (k  =  1) 


2  32v  .  3v  A  ,2  3v. 

r  -T  t  2r  r-  0  r  -  (r  —) 


This  is  satisfied  by  the  equation 
v  =  w  -  c[l  -  (a/r)]  where  w  and  c  are  constants. 
E.  EQUATION  SOLUTION 

rL  ♦  (It  k)r  ~  -  (1  -  k)v  =  0 
3r2  3r 


Lira  [v]  =  w  ,  Lim  [v]  *  0 

r  -*  a  r  •*  • 


Use  of  the  substitution  r  *  aea  gives 

.2  ■> 

3  v  a  .  }v  .  . 

— =•  ♦  k  —  -  (1  -  k)v  ■  0 

*a2  * 


Lim  [v]  »  w  ,  Lim  [v]  «  0  , 

a  -*  0  7  -*  • 


which  has  the  standard-fora  solution 


v  «  Ae"**  ♦  Bo' 


(k  -  1)« 


A  ♦  B  ■  w  ,  k  >  1  , 

which  in  terns  of  the  original  radial  coordinate  r  becomes 

.  .a.  .  .a.k— 1 

v  •  A  •  (-)  •  B  •  (-) 

r  r 

A  ♦  B  «  w  ,  k  >  1 
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F.  RAOTATIONAL  KINETIC  ENERGY 

Let  h  be  the  vertical  thickness  of  a  disk-shaped  slice  made 
through  the  tornado  parallel  to  the  ground,  and  let  E  be  the  ro¬ 
tational  kinetic  energy  of  this  slice  extending  radially  from 
r  =  a  to  r  =  -. 


/Ot  n  V  #00  2 

k-t-  •  2urh  •  dr  =  irph  j  rv  dr 


=  -ph  i  r  •  [-j  »  a  rr3  ir 

a  P  r  r 

k  2  2k  -  2 

,  r.  2,  .  2ABa*  _2  -  k  .  Ba*  4  -  2k,- 

=  nph  [Aa  Ln  r  +  2'  -  %  r  ♦  "  4 "IT" 2 -  r  -*a 


It  is  apparent  from  the  result  that  E  will  be  finite  if  and  only 


if  A  *  0  and  k  >  2.  Consequently, 


2  2  2  2 
wphB  a  ira  hpw 


ira  how 

STk  -  5) 


k  >  2  ,  V  a  W  •  (|-)k  "  1 


G.  ANGULAR  MOMENTUM 


J  *  prv2*rh*dr  *  2wph  r^vdr 


2«eh  £ 


k  -  1  #  .  1  -  k  . 

1  •  r  •  r  •  or 


k  - 1  IV  - k“  * 

•  2»oh  w*  •  ^ 

“  -  a 

It  is  apparent  from  this  result  that  J  will  be  finite  if  and 
only  if  k  >  4.  Consequently, 

.  .  2«phu3w  i_  .  .. 

J  *  —  »  k  >  4  . 
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H.  GRAPHICAL  RESULTS 

The  following  page  (Figure  21}  shows  in  graphical  form  the 
circling  wind  speed  outside  a  tornado  funnel,  calculated  from  the 
tangential  velocity  formula  for  k  «  4,  k  s  6,  and  k  =  8.  Since  the 
ai  ilar  momentum  of  the  tornado  is  finite  only  if  k  >  4,  it  is  seen 
that  in  all  admissible  cases  the  circling  wind  speed  falls  to  less 
than  ten  percent  of  its  funnel-wall  value  beyond  a  distance  of  three 
funnel  radii. 
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.1  .  CONCLUSIONS 

Tliu  x.jot:  mathematical  equations  for*  the  three-dimensional 
viscous  and  compressible  air  flow  in  the  proximity  of  a  tornado 
funnel  are  extremely  complex,  and  require  a  large  investment  of 
further  analytical  and  computer  effort  in  order  to  develop  fully 
adequate  solutions.  The  preliminary  analysis  presented  here  deals 
only  with  two-dimensional  laminar  and  incompressible  flow  in  a 
thick  disk  aligned  co-axi'lly  with  the  funnel  axis. 

A  rather  peculiar  formula  for  the  tangential  air  velocity 
resulted  from  the  solution  of  the  equations  of  the  restricted  system. 

The  circling  wind  speed  was  shown  to  decrease  inversely  with  the 
(k  -  1)  power  of  the  distance  from  the  funnel  axis,  where  k  is  a 
rather  unique  dimensionless  constant.  The  denominator  of  the  d;- 
mensionless  ratio  k  is  equal  to  the  kinematic  viscosity  u/p  of  the 
air,  and  the  numerator  is  equal  to  the  algebraic  product  of  the  funnel 
radius  and  the  radial  speed  of  air  flow  into  the  funnel.  It  was 
found  that  rotational  kinetic  energy  of  the  soinning  air  disk  cannot 
be  finite  unless  k  >  2,  and  that  the  angular  momentum  of  the  spinning 
diks  cannot  be  finite  unless  k  >  4.  The  constraint  that  k  >  '* 
requires  that  the  circling  wind  speed  must  decrease  very  rapidly 
with  distance  from  the  tornado  funnel  well. 

Inspection  of  photographs  made  of  actual  tornados  indicates  that 
the  assumption  of  laminar  flow  is  probably  not  a  serious  restriction 
except  in  the  small  region  where  the  tornado  funnel  touches  the 
ground.  The  assumption  of  incompressibility  is  more  serious  in  the 
vicinity  of  the  tornado  wall,  where  the  circling  wind  speed  Is 
'roquontly  said  to  approach  half  the  speed  of  sound. 

The  more  elaborate  study  should  concentrate  first  on  the  removal  of 


the  incompressibility  restriction  and  then  take  into  account  the 
vertical  component  of  air  flow  as  influenced  by  gravity  and  the 
equation  of  continuity. 
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